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When the rare earth mononitrides (RENs) first burst onto the scientific scene in the middle of last 
century, there were feverish dreams that their strong magnetic moment would afford a wide range of 
applications. For decades research was frustrated by poor stoichiometry and the ready reaction of 
the materials in ambient conditions, and only recently have these impediments finally been overcome 
by advances in thin film fabrication with ultra-high vacuum based growth technology. Currently, the 
field of research into the RENs is growing rapidly, motivated by the materials demands of proposed 
electronic and spintronic devices. Both semiconducting and ferromagnetic properties have been 
established in some of the RENs which thus attract interest for the potential to exploit the spin 
of charge carriers in semiconductor technologies for both fundamental and applied science. In this 
review, we take stock of where progress has occurred within the last decade in both theoretical 
and experimental fields, and which has led to the point where a proof-of-concept spintronic device 
based on RENs has already been demonstrated. The article is organized into three major parts. 
First, we describe the epitaxial growth of REN thin films and their structural properties, with an 
emphasis on their prospective spintronic applications. Then, we conduct a critical review of the 
different advanced theoretical calculations utilised to determine both the electronic structure and 
the origins of the magnetism in these compounds. The rest of the review is devoted to the recent 
experimental results on optical, electrical and magnetic properties and their relation to current 
theoretical descriptions. These results are discussed particularly with regard to the controversy 
about the exact nature of the magnetic state and conduction processes in the RENs. 



I. INTRODUCTION 

The rare-earth mononitrides were first investigated in 
the 1960s, when technological developments overcame 
the problems faced in separating the chemically similar 
members of the lanthanide series. That early wo rk was 
discussed in reviews from the 1960s and 1970sP2l i n the 
following period, from 1980 to 2005, there was a low level 
of interest, with further reviews appearing in the 1990s 
and in the early years of this centuryi^nSI l n the past few 
years the literature on the subject has grown rapidly, 
based on a combination of breakthroughs in theoreti- 
cal developments and the facility to grow stoichiomet- 
ric epitaxial films. The rare-earth nitrides (RENs) show 
promise in applications as diverse as spintronics, infrared 
(IR) detectors and as contacts to III-V compounds, with 
now about a dozen laboratories worldwide reporting the 
growth and study of REN thin films. 

The rare-earths, with atomic numbers from 57 (La) 
to 71 (Lu), comprise the elements across which the 
4/ orbitals are filled. They have atomic configurations 
LYe]6s 2 4/™ , with n varying from 1 for Ce (0 for La) to 
14 for Lu, and certain of the series have an additional 



5d electron (La, Ce, Gd, Lu). Their most common ionic 
charge state is 3+, with the 4/ levels spanning the Fermi 
energy. They are the only stable elements with more 
than marginally filled /-shell electronic orbitals, and as 
a consequence they are the elements with the largest spin 
and orbital moments. In ordered solids they contribute 
to the most strongly ferromagnetic materials, a contri- 
bution that has ensured their utility in technologies that 
require strong permanent magnets. Despite their name 
they are by no means rare, with the exception of prome- 
thium, which has no stable nuclear isotope. They are 
found in the Earth's crust at concentrations exceeding 
that of Ag, Cd and Hg, and similar to Ge and AsP 

The magnetic states of most rare-earth monopnictidcs 
(RE-V) have been known for some decades, after stud- 
ies undertaken already in the 1960s. The heavier pnic- 
tides were found to be antiferromagnetic, but in contrast 
the nitrides are almost all ferromagnetic; the magnetic 
properties are reported in previous reviews The RE- 
V adopt magnetic order only at cryogenic temperatures, 
with the highest Curie temperature (Tc), for GdN, of 70 
KpEl The current interest in the potential of spintronic 
devices has raised the level of urgency in the exploration 
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of intrinsic ferromagnetic semiconductors, of which the 
RENs offer a rich set of examples. Yet their potential has 
been poorly explored until recently, despite EuO having 
become a well-known intrinsic magnetic semiconductorP 
Unlike the dilute magnetic semiconductors (DMS) the 
RENs do not rely on the presence of foreign ions, nor 
on a huge hole concentration that prevents independent 
doping control. Thus in principle this permits controlled 
doping independent of the ferromagnetism opening new 
possibilities for spintronic devices and fundamental spin- 
transport research. However even within the DMS sce- 
nario the RENs hold some promise, with reports that 
Gd leads to room-temperature ferromagnetism with an 
unexpectedly large moment when it is introduced as a 
dilute impurity in GaNP^l A thorough understanding of 
the corresponding GdN is a prerequisite to understand 
this proposed DMS system. 

In contrast with the magnetic properties, the transport 
properties and electronic band structures of the RENs 
were until recently much less certain. Even for ostensibly 
identical compounds one can find reports of the conduc- 
tivity ranging from insulating to metallic. The source of 
the uncertainty lies in two issues that impact strongly 
on their achievable stoichiometry: a propensity for rapid 
oxidation when exposed to air and for the formation of 
N vacancies (Vjv). The former requires that any thin 
film must be protected by a passivation capping layer, 
while the latter is related to the small formation energy 
of Vm that ensures these will be present at the level of at 
least 1% in any film grown substantially above ambient 
temperature. 

Among the earliest of the recent theoretical discussions 
there was a prediction that the RENs would display a se- 
ries of contrasting conductive states, ranging from semi- 
conductors to semimetals, and including half metalsP^ 
This range of interrelated magnetic and conducting prop- 
erties is then an obvious testing ground for prototypical 
spintronics structures. Following that prediction there 
have been a number of laboratories that have initiated 
programmes to grow polycrystalline and epitaxial films, 
so there is the hope that some spintronic structures can 
ultimately be explored. That hope is still some way off, 
waiting for improved techniques for the growth and cap- 
ping of thin films of these severely reactive materials, 
although already now there is a report of a GdN-based 
spin filter 

Against this background there has been for some time a 
theoretical interest in the rare-earth monopnictides. The 
localized atomic-like properties of the open-shell 4/ elec- 
trons defy the standard density functional band structure 
theory approach. Because the RENs form a family of 
materials with the same simple rocksalt structure, they 
have become a useful testing ground for new theoreti- 
cal developments dealing with the strongly correlated 4/ 
electrons. That work is impeded by the paucity of re- 
liable experimental data, especially about the electronic 
state of the nitrides; there is relatively little too inform 
the theoretical work. The recent rapid expansion of ex- 



perimental studies, especially on well-ordered films, has 
begun to provide data with which to tune the treatment 
of strong correlation. 

In this review we will focus on the recent developments 
on the rare-earth (RE) nitrides. However, occasionally, 
we will also include closely related materials, such as 
other RE-V, as well as some RE chalcogenides and ox- 
ides, such as EuO. For example, because the RE shell 
is half-filled in both EuO and GdN, it is of interest to 
compare them. We should mention that the theoreti- 
cal properties of RE-V were reviewed not so long ago 
by Duan et alP Nonetheless, the present review, focused 
on RENs will be complementary, particularly because of 
their contrasting magnetic behaviour the recent experi- 
mental studies based on advances in epitaxial thin-film 
growth. 

We have prepared the present review with an eye to 
the potential exploitation of this class of material in var- 
ious applications. Their potential depends on specific 
aspects of their crystal structures and electronic prop- 
erties, and although these will be discussed in detail in 
later sections, in the following section we give a brief 
description of the most obvious applications. Specific 
attention will be paid to issues such as their epitaxial 
compatibility. Following that we will cover recent ad- 
vances of, in Section III epitaxial thin-film growth and 



in Section [TV] theoretical descriptions. The magnetic and 
electronic/magnetoelectronic advances will be discussed 
in Sections [V] and VI respectively, and Section VII is a 
summary. 



II. THE BASIS OF POTENTIAL 
EXPLOITATION 

The RENs form in the face-centered (FCC) cubic NaCl 
structure with lattice constants ranging from 5.305 A for 
LaN to 4.76 A for LuN, in total a 5% difference across the 
series and less than 0.5% between nitrides of neighbour- 
ing atomic species. There is clearly potential for epitax- 
ial growth of custom-designed heterostructures, including 
superlattices, and even for controlled strains to be intro- 
duced. In later sections we will discuss the complemen- 
tary electronic properties and strongly contrasting mag- 
netic behaviours of the RENs, traits which immediately 
suggest them as the basis for a variety of spin-dependent 
devices. It is essential for any such devices that tech- 
niques are developed for the growth of well-ordered epi- 
taxial structures and that a thorough understanding is 
reached of the electronic band structures and magnetic 
behaviours of the RENs, both in bulk and in thin epitax- 
ial layer forms. 

The strong exchange interaction results in a significant 
spin splitting of a few hundred meV in both band edges, 
with the majority spin having the lower energy in the 
conduction band and the higher energy in the valence 
band. Thus carriers at both edges, electrons and holes, 
are in majority-spin bands; the minority spin edges are 



3 



unoccupied at ambient temperature. Any device, such as 
a diode, transistor or filter that requires doping or accu- 
mulating carriers into the band edges will involve trans- 
port of carriers with only majority spin state. Based on 
this strong exchange splitting, a REN-based spin-filter 
Josephson junction has been recently demonstrated.^ A 
thin layer (5 nm) of GdN acting as a spin-dependent 
tunneling barrier is placed between two NbN supercon- 
ductor contacts, with the resulting spin-filter efficiency 
estimated to be about 75% at 4.2 K. With a strong mag- 
netoresistance over a broad field range in its ferromag- 
netic state, as well as a small coercive field, GdN is an 
appropriate material to explore for use in magnetic field 
sensors. The recent enhancement of the Curie tempera- 
ture up to 200 K in N-deficient GdN films would allow 
operation of the sensor well above the boiling point of 
liquid nitrogen! 13 * 14 ^ 

With magnetic states that vary strongly across the se- 
ries and coercive fields depending strongly on the growth 
conditions, the RENs are of specific interest for non- 
volatile magnetic memory elements such as MRAMs. In 
particular, SmN is the only known near-zero-moment fer- 
romagnetic semiconductor, with an enormous coercive 
fieldp^ and, combined with GdN which has a coercive 
field some three orders of magnitude smaller, they may 
form an ideal hard- and soft-ferromagnetic pair with po- 
tential in memory elements. 

The ongoing effort for the epitaxial growth of RENs is 
further driven by the potential to integrate REN-based 
devices with group-Ill nitrides to develop new function- 
alities combining both families. The narrow band gaps 
of the RENs, with optical absorption edges lying near 1 
eV and absolute gaps on the order of one half of that, are 
interesting for IR detectors. Thus the properties of the 
RENs are complementary with those of the wide band- 
gap group-Ill nitrides, and a heterojunction involving the 
two semiconductors might have very attractive properties 
for multi-wavelength photonic devices. In addition, as 
spin-polarised carrier injection cannot be accomplished 
efficiently from metals into a semiconductor, a GdN layer 
could be regarded as a spin injector in GaN-based tran- 
sistors or diodes. 

Although the Si-RE reaction is so rapid that it prevents 
the epitaxial growth of RENs directly onto Si, they can 
be grown as polycrystalline films at ambient tempera- 
ture, which may provide useful injection of spin selected 
electrons. The use of A1N and GaN buffer layers will 
be seen to permit epitaxial REN/Si integration; other 
buffers have yet to be explored.^ 

Very recently there has been an ab initio-b&sed sug- 
gestion that GdN will have a ferroelectric ground state 
under 3% compressive in-plane strain, and a similar re- 
sult can be expected for the other RENs.^ The lattice 
constant variation across the REN series permits a thor- 
ough investigation of the possibility. If the prediction 
is verified they would form interesting electromagnetic 
multifcrroics. 

Thus among the most interesting aspects of the series is 



their epitaxial compatibility coupled with their contrast- 
ing magnetic and complimentary electronic properties, 
showing promise for a wide range of possible spintronic 
and electronic structures. It is then of special interest 
to provide a much fuller description of the magnetic and 
electronic properties of the entire series. To date it is 
only GdN that has been subjected to very thorough ex- 
perimental investigation, as will be realised very quickly 
on reading this review. 

Finally, there are reports of catalytic^ and large 
magneto-caloric^HUl effects in the RENs, suggesting 
them as promising candidates for magnetic refrigeration. 
Thin films are unlikely to contribute to these technolo- 
gies. 



III. CRYSTAL STRUCTURE AND EPITAXIAL 
GROWTH OF RARE-EARTH NITRIDES 

One of the major hurdles in growing REN thin films, 
either epitaxial or in polycrystalline form, is their propen- 
sity to form V jy an d to decompose in air into RE ox- 
ide/hydroxides This is the reason for some of the 
long standing controversies concerning their electronic 
structure, transport properties and magnetic properties. 
The growth of RENs in both thin film and bulk form 
dates back more than seventy years and quite a number of 
processes have been tested, but none was found fully sat- 
isfactory at that timeP It is only recently that high qual- 
ity epitaxial thin films have been achieved, mainly thanks 
to the continuous development of film growth techniques, 
such as ultrahigh vacuum (UHV)-based methods. In this 
section we describe the recent developments in the field 
of epitaxial growth of REN thin films. The growth of 
bulk materials will not be treated in this review. Poly- 
crystalline thin films will be mentioned briefly, primarily 
because their lower growth temperature offers some con- 
trol over the formation of VjvP 



A. Crystal structure 

The RENs adopt a FCC NaCl structure (space group 
Fm-3m (225)) as shown in Figure [I] and in all NaCl- 
type RENs with the exception of CeN (tetravalent) the 
cation is trivalent. Each RE atom is coordinated by 6 
nitrogen atoms and, conversely, each nitrogen atom is 
coordinated by 6 RE atoms. The large difference in elec- 
tronegativity between nitrogen (3.0) and RE (1.1 to 1.5) 
leads to a strong affinity and a predominantly ionic char- 
acter (more than 50%) of the RE-N bonds. 28 High quality 
single crystals have proven to be exceedingly difficult to 
prepare and as mentioned above there remains much con- 
troversy about various experimental results, and many 
fundamental material parameters have still not been re- 
liably measured. 

Table |T] shows some of the crystal lattice parameters 
for the RENs available in the literature, which range 
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FIG. 1. The common rocksalt crystal structure of the RENs. 
The large spheres represent the RE cations while the small 
spheres are the N anions. 

from about 5.305 A to 4.76 A for LaN to LuN, re- 
spectively. The lattice constants of the RENs decrease 
with increasing 4/ occupancy as expected from the de- 
crease of cation sizes across the series. The theoretical 
lattice constants for the entire REN series determined 
from first-principles calculations are shown in Tabic [I] 
and differ from most of the experimental lattice constants 
by about 1-2% The bulk moduli calculated following 
the same computational approaches are also displayed in 
Table [I] with in this case larger discrepancy among the 
calculation methods. However, the calculated values are 
still comparable to experimental values determined for 
the 5/ nitrides UN (203GPa) and ThN (175GPa)p2l and 
for the group Ill-nitrides A1N (209GPa), GaN (183GPa) 
and InN (133GPa)P^ No experimental data for the ther- 
mal expansion, thermal conductivity or elastic stiffness 
coefficient, including Young's modulus, have yet been ob- 
tained, with the exception of CeN. For the latter, the 
bulk modulus of 153 GPa has been determined by X- 
ray diffractiorfSl while the hardness and elastic modu- 
lus of epitaxial CeN(OOl) films were determined from 
nanoindentation measurements to be 15.0 ±0.9 GPa and 
330 ± 16 GPa, respectively.^ A full computational study 
of the elastic properties of CeN has been performed by 
Kanchana et alP2 and the results are in close agreement 
with the experimental values. More recently, several 
groups have theoretically explored the elastic properties 
and hardnesses of PrN™ TbJj^H and HoNp3 and Yang 
et al. have studied the entire REN series^ and compared 
their results with the literature. Furthermore, linear di- 
electric constants of the RENs derived from the dielectric 
theory of chemical bonds for solids are also displayed in 
Table EF 



TABLE I. Experimental lattice constants, calculated lattice 
constants, bulk moduli (B) and dielectric constants (e) of 
rare-earth nitrides. 



REN 


Expt. 




Calc. 






a(A) 


a(A) 29 


B{GPa ^mmM 


,110 


LaN 


5.30fP 


5.38 


130, 122 


7.07 


CeN 


5.022^1 


4.90 


210, 121 


7.01 


PrN 


5.13EP 


5.29 


140, 121, 129.14 


6.96 


NdN 


5.132^1 


5.24 


140, 197 


6.90 


PmN 




5.19 


150, 114 


6.85 


SmN 


5.03P 


5.10 


180, 127 


6.79 


EuN 


5.01<E11 


5.14 


110, 114 


6.74 


GdN 


4.97#H 


5.08 


150, 111 


6.68 


TbN 


4.92P 


5.05 


150, 241, 155.53 


6.63 


DyN 


4.89P 


5.03 


160, 121 


6.57 


HoN 


4.86P 


4.98 


170, 138, 137.9 


6.51 


ErN 


4.845P 


5.00 


160, 56 


6.46 


TmN 


4.8P 


4.90 


190, 138 


6.40 


YbN 


4.781 1331 


4.79 


190, 136 


6.35 


LuN 


4.7<P 


4.87 


170, 183 


6.29 



B. Epitaxial growth of REN 

The progress in studies of the RENs over the last 
few years was achieved using UHV- based methods 
such as molecular bea m epitaxy (MBE), 16 46 ~ 5 1 pulsed- 
laser depo sition (PLD) j^ 50 ' 52 ' 53 1 and dc/rf magnetron 
sputtering! 54 ^ These are the techniques of choice for de- 
position of RE-based materials because the high- vacuum 
and inert-gas environment helps to ensure material purity 
and interface quality and these methods are subsequently 
the most common used in this field. Among these UHV- 
based techniques, the recent advances in the growth of 
GdN and EuN by MBE tend to show that this is probably 
the best growth method to ac hieve rapid improvement in 
the quality of REN thin films! 16 ' 47 1 49 ^ 1 It is worth men- 
tioning that the purity of the as-received RE charges or 
targets is far from the one commonly used for the growth 
of conventional III-V semiconductors; typically the pu- 
rity of the RE source is about 99.99% by weight with the 
main impurities being oxygen, carbon, nitrogen, calcium 
and iron. 

While solid sources arc commonly used for the RE ele- 
ments, a wide range of options are available for nitrogen 
sources/precursors. The most straightforward is to use 
pure N2 gas by taking advantage of the catalytic break- 
down effect of the RE atoms o n the molecular nitrogen. 
So far epitaxial growths of Gdl>J 16 1 50 1 and SmP^ as we n 
as polycrystalline films of DyN, ErN, LulNP^ have been 
demonstrated. The growth of EuN must be achieved un- 
der the presence of radical nitrogen, for example using 
low energy nitr ogen ions from a Kaufmann ion source 
or N 2 plasmaP^ l 50 * 52 ! The other nitrogen precursors used 
for the growth of REN are those commonly used for the 
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growth of GaN: ammonia (NHsj^, which is decomposed 
on the surface of the substrate by pyrolysis, and nitrogen 
plasma obtained either by r adio- frequency or electron cy- 
clotron resonance! 47 ' 48 ' 52 ' 5 ^ Epitaxial growth of the RENs 
is still at too early a stage to conclude on the best nitro- 
gen source in terms of material quality. 

The growth conditions of CeN, GdN and SmN un- 
der NH3, pure N2 and low energy nitrogen ions are in 
fact close to those currently used for the III-V nitrides. 
They are grown with an excess of nitrogen species with 
respect to the RE flux to avoid Vjv formation and/or 
metallic RE clusters. Nitrogen ions to Ce ratio of 
and pure or NHj^Il to Gd/Sm ratios larger than 

100 have been reported. Using these conditions, stoi- 
chiometric films are achieved and the material properties 
are of high quality. Epitaxial growth of GdN has also 
been performed usin g an N2 plasma cell with a slight 
exce ss of Gd by MBlP^ and in a large excess of N2 by 
PLDpJESl The growth of EuN is somewhat intriguing as 
epitaxial films have been obtained when the growth oc- 
curs in an adsorption-controlled growth regime, where 
the Eu flux is set to be much higher than the N 2 flux and 
a high substrate temperature (~800°C) is maintained to 
re-evaporate any excess of EuP^EuN has also been grown 
by PLD but the authors do not report on the RE:N ratio 
while they point out that high growth temperatures, up 
to 860°C, are required to achieve epitaxial films An op- 
timised growth temperature window around 600±50°C 
was found for SmN, 50 again in a range of growth tem- 
peratures and pressures where any excess of Sm can be 
re-evaporated from the substrate. 

While it seems there is a consensus about the growth 
temperature for EuN and SmN thin films, the situation 
is far from clear about the one to use for GdN. Gerlach 
et alP^ reported a growth temperature of 750° C, com- 
parable with those used later on by Natali et alP^ for 
the growth of GdN under pure N 2 (650-750° C) or by 
Ludbrook et al. (700-850° C) by PLD using a N 2 plasma 
source.^ No clear significant temperature dependence of 
the structural quality of the films was found between 700 
and 850°C. Temperatures much lower, such as 500°C, 
have been used for the growth of epitaxial GdN films by 
reactive radio frequency magnetron sputtering! 5 - 5 -' GdN 
layers have also been grown at 450° C by MBE using an 
N 2 plasma cell and NH 3 as the nitrogen precursor.^ The 
authors mentioned that deposition at higher tempera- 
tures did not yield GdN films. It is worth mentioning 
that Gd has a vapor pressure of 1.3 10 -4 mbar at 1175°C, 
and, therefore, the growth process does not involve the 
re-evaporation of the excess Gd. Only one paper reports 
on the epitaxial growth of CeN, for which the growth 
temperature was 700° CP^ 

Although metal-organic chemical vapor deposition 
(MOCVD) is the technique of choice for the growth of 
group-Ill nitrides, for the RENs the weak point of the 
technique has been for a long time related to the ab- 
sence of an efficient RE precursor. Only recently have 
GdN and DyN polycrystalline films been achieved using 



guanidinato-complexes o f Gd a nd Dy as precursors! 57 ' 58 ' 
Physical vapor depositio n 1 59 ' 60 ' and more recently chemi- 
cal vapor depositiorPS as well as plasma-enhanced atomic 
layer depositiorP^l have been used to grow polycrystalline 
GdN films. 



C. Substrates and capping layers 

One of the major difficulties which has hindered REN 
epitaxial growth is the lack of native substrates. A 
plethora of suitable material that is more or less matched 
with the RENs has been employed. Below we give a brief 
description of each of these substrates. 

Most research in the past has selected (100) oriented 
substrate surfaces for REN epitaxy. This is supported 
by the fact that the RENs adopt a FCC (NaCl) struc- 
ture. Historically the first epitaxial growths of RENs, 
CeN by T.-Y. Lee et alPin 2003 and then GdN by Ger- 
lach et alPin 2007, have been performed on MgO(100) 
substrates with a lattice mismatch of about +18.7% and 
+19.2% for GdN and CeN, respectively (Figure^. The 
film/substrate epitaxial relationships were demonstrated 
to be (001)REN||(001)MgO and [100]REN|| [100] MgO. In 
both cases the films were shown to be of high crystalline 
quality, but in the case of GdN little was reported on 
the magnetic or transport properties and it has been ob- 
served that films thicker than 60 nm have very rough 
surfaces. The most closely lattice-matched substrate for 
REN epitaxy would be YSZ (100) with a lattice param- 
eter of 5.125 A, nearly matching the lighter RENs as 
show n in F igure [2] Epitaxial growth on YSZ substrates of 
GdN j5om] SmNp an d EulsP using PLD and then later 
EuN by MBlpSlhave been reported. It is worth mention- 
ing that an oxide layer, RE 2 03, is formed at the interface 
between the substrate and the REN films, likely due to 
the strong affinity of the RE for oxygen and the mobil- 
ity of oxygen in YSZp^ For economic and technological 
reasons it would be a real advantage to demonstrate the 
growth of REN films on silicon substrates. The REN lat- 
tice constants are within about 10% of that of Si (Figure 
[2J, but, as mentioned previously, silicide formation at the 
Si/REN interface is still a major issue to overcome and 
has so far prevented epitaxial growth. 

The possibility to take advantage of the hexagonal 
symmetry of the {111} plane of the REN rocksalt struc- 
ture to investigate growth on c-plane (0001) wurtzite 
GaN surfaces has been propo sed b y Scarpulla et al.P 
and then later on A1N surfaces! 16 ' 55 ' Earlier work by Shi- 
momoto et al. reported the growth of EuN as a buffer 
layer on MgO(lll) and AL2O3 substrates for the subse- 
quent growth of InNP^l Figure [2] shows the lattice mis- 
match between (111) REN lattice and c-plane (0001) 
wurtzite A1N, GaN and InN lattices as a function of the 
RE atomic numberP^ This graph shows that RENs can 
be grown nearly lattice-matched on an InN substrate, 
but its lack of commercial availability make GaN and 
A1N more attractive for the moment. The 7.5-15% lat- 
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in the past, including metal lic layers such as W,^l Q r JSH 
Cu^ Tar P 1 and NblN P ^I and i nsulator films such as 
YSZ jHl GaNj^ S | 47 | B4 | b5 1 Al^PlM and MgF 2 PPolycrys- 
talline A1N and GaN are the most commonly used cap- 
ping layers which can be attributed to their transparency 
allowing optical measurements, their ease of growth, and 
good chemical stability over time. 



D. Structural Properties 
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FIG. 2. Calculated lattice mismatch between RENs and sub- 
strate materials. Top: (001) REN and MgO(OOl), YSZ (001) 
and Si(001). Bottom: (111) REN c-plane (0001) wurtzite 
AIN, GaN and InN lattices as a function the RE atomic num- 
ber. From Reference 1501 



tice mismatch between the RENs and either AIN or GaN, 
though relatively severe, is comparable to that in het- 
eroepitaxial systems.^ The growth of the REN (111) 
plane on a (0001) surface resulted in two rotational vari- 
ants of the grains due to the lower 3-fold symmetry of 
GdN compared to the 6-fold symmetry of the group- 
Ill nitride surface.^ The in-plane epitaxial relationship 
for the two rotational variants are [l-10]GdN|| [10-10]GaN 
and [-1010]GaN, in addition to (lll)GdN|| (0001)GaNP 
The quality of the films of Refp3 was high enough to ad- 
dress the properties of the RENs with relative confidence. 
To date most of the laboratories have selected group- 
Ill nitride substrates for REN epitaxial growth. This 
choice is supported by the fact that these substrates are 
widely available with a good crystal quality, due to their 
transparent nature, and ease of handling and pre-growth 
cleaning. Their use is also driven by the possibility to 
develop new functionalities combining nitride semicon- 
ductors and REN materials. Additionally, the possibility 
to grow GdN on a 100-nm thick AIN buffer layer on top 
of a silicon substrate has been demonstrated^ and may 
pave the way for the integration of RENs into mainstream 
silicon technology. 

Due to their decomposition in air, REN films must be 
passivated with an effective capping layer to avoid re- 
action with the atmosphere. A series of polycrystalline 
or amorphous capping layers have been tried successfully 



In spite of efforts to improve crystal quality, only little 
is known about the structural defects in the RENs. This 
is mainly due to their instability in air which makes high 
resolution transmission electron microscopy challenging. 
However, for some of the RENs, such as GdN, the oxida- 
tion rate is sufficiently s low to allow cross-section scan- 
ning electron microscopy! 47 * 55 * 6 -^ As most of the RENs are 
prepared by heteroepitaxy on III-V nitride substrates, 
the most prevalent structural defects should be thread- 
ing dislocations. In addition, a misfit type of disloca- 
tions lying at the interface to relieve the misfit strain in 
the large lattice mismatched growths is expected what- 
ever the growth orientation, (001) or (111). Reflection 
high energy electron diffraction (RHEED) has been used 
for measuring the plastic relaxation during the growth of 
GdN on AINpllt has been observed that after ~2.5MLs 
the lattice constant of GdN starts to increase, while the 
RHEED pattern remains streaky, and a fully relaxed 
GdN layer is obtained after only 6MLs. RHEED has 
also been useful to monitor the presence of RE2O3 at the 
interface during growth on a YSZ substrate. 50 ^Interest- 
ingly it has also signalled the presence of a (2x2) surface 
reconstruction after the growth of (001) EuNp^ while 
no surface reconstruction was observed for (lll)-oriented 
films. In the case of growth onto group Ill-nitride sur- 
faces, planar defects in the form of twin boundaries ro- 
tated by 60 degrees with respect to each other are present 
in the films These are separated by grain bound- 
aries and/or antiphase domains, meaning that single- 
crystalline films will only be possible when one variant 
can be suppressed. 

The traditional figure of merit used to assess the qual- 
ity of the films is the X-ray diffraction (XRD) coherence 
length determined from 8-26 and rocking curve scans. 
Full-width half maximums (FWHMs) for the symmet- 
ric (002) rocking curves of 2.5° and 1.4° have been re- 
ported for (OOl)-oriented GdN grown on YSZ(001^and 
MgO(00ipl respectively, while a larger FWHM, 2.16°, 
for the symmetric (002) rocking curve of (OOl)-oriented 
CeN grown on MgO (001}P^ has been measured. Nar- 
rower linewidths are found for films grown in the (111) di- 
rection. Scarpulla et alP^ reported values with a FWHM 
for the symmetric (111) peak of 0.251° and 0.321° for 
GdN grown on GaN using an N2 plasma or with NH3 
as a nitrogen source. Such values lead to in-plane X-ray 
coherence lengths up to 100 nm. While the structural 
quality of these films is higher than for any other re- 
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ported films, further optimisation of the growth parame- 
ters is still required to reduce the degree of twinning and 
mosaic spread. 



IV. RECENT THEORETICAL ADVANCES 

The theoretical discussion starts with a review of some 
basic textbook notions about open-shell 4/ systems, such 
as Hund's rules, in Sec. IV A| Next, we will briefly review 
the shortcomings of standard density functional theory 
for 4/ systems. In Sec. |IVB| we will review the early band 
structure work on RE-V, which used the simplest possible 
approach, namely, treating the 4/ electrons as a partially 
filled core state. Then we will turn to each of the main 
new theoretical approaches for dealing with the strongly 
correlated 4/ electrons: the self-interaction correction 



approach (Sec. |IVC| , the LSD A+E/ method (Sec. |IVD| , 

GW theory (Sec. 



IVE) 



the hybrid functionals (Sec. 
IV F ) , and the dynamical mean field theory (DMFT) ap- 
proach (Sec. IV G). All of this will mainly be focused 



on the electronic structure. One of the main questions 
to be answered is whether the RENs are semimetals or 
semiconductors. 

A large fraction of the existing theory work is devoted 
to understanding the magnetic properties, in particular 
why GdN is ferromagnetic while the other pnictides are 
antiferromagnetic. On the other RENs, there are far 
fewer detailed studies of the magnetic properties from 
a band structure point of view, although there is a large 
body of older work dealing with the complex magnetic 
ordering phenomena in these materials. Besides the mag- 
netism there is also some interesting theory work on the 
optical and lattice dynamical properties of the RENs. 
Therefore, we discuss the magnetism (Sec. [v]), electronic 
and optical properties (Sec. VI) and vibrational proper- 



ties (Sec. VIC I in separate sections. 



A. Basic properties of 4/ electrons. 

The atomic-like properties of 4/-electrons arise from 
their localized nature. The high angular momentum 
I = 3 of an / electron results in a strong centrifugal 
barrier 1(1 + l)/r 2 which keeps the electrons separated 
from the other valence shells, like the 5d and 6s elec- 
trons. Because the 4/ radial wave function has no node 
it is rather localized. However, it increases as r 3 near the 
origin and is thus peaked away from the nucleus, but at 
the same time at a much smaller radius than the outer 
6s and 5c? electrons which need to stay orthogonal to the 
lower s and d shells. In contrast, when we consider the 
actinides, the 5/ shell is more spread out and behaves 
less atomic-like. Of course, how localized or band-like 
the electrons are depends on the crystal structure and 
on the separation of the atoms. Hence, there is some in- 
terest in early lanthanides, whether under pressure they 



can convert from an atomic-like to a band-like behavior 
of the 4/ electrons. 

The atomic localized character implies a strong 
Coulomb interaction between the electrons as we grad- 
ually fill the 4/ shell. The basic atomic physics is 
thus dominated by the electron-electron interaction. In 
atomic physics, the basic question to answer is then, what 
are the quantum numbers that determine the ground 
state of the multi-electron system consisting of an f N 
system? 

This question has been studied long ago in terms of 
the different coupling schemes: L • S coupling, J • J cou- 
pling, etc., and is for example summarized in the well- 
known book by Condon and Shortley.S2 If the spin-orbit 
coupling is smaller than the electron-electron interaction, 
one expects the L • S coupling scheme to be valid. The 
many-electron wave function is a linear combination of 
Slater determinants with a given total orbital and spin 
angular momentum L and S and these then couple to 
give a total angular momentum J. In that case, the good 
quantum numbers are L, S, J, Mj. Here L, S, J define the 
eigenvalues of L(L + 1)H 2 , etc., of the angular momentum 
operators L 2 , S 2 , J 2 . These determine the magnitude of 
the corresponding angular momenta. The Mj determine 
the eigenvalues of the total angular momentum along one 
chosen quantization axis, J z . As is well known, quantum 
mechanics allows us only to measure the magnitude and 
one component of the angular momentum, but not the 
various components simultaneously. Determining all pos- 
sible combinations of these quantum numbers for a given 
number of / electrons determines the so-called multiplct 
splitting terms. 

The above theory for free atoms results in the well- 
known Hund's rules, which determine which of the an- 
gular momentum quantum numbers result in the lowest 
energy. These rules are: (1) first the total spin S should 
be maximized, (2) then, for a given maximal spin, the to- 
tal orbital angular momentum L should be maximized, 
(3) if the shell is less than half filled, then the lowest 
energy is found for J — \L — S\, while for a more than 
half-filled shell, J = L + S. In other words, the third 
rule results from the spin-orbit coupling and says that 
the orbital and spin momenta oppose each other for less 
than half-filling and are in the same direction for more 
than half-filling. As an example, take the Eu 3+ ion with 
6 / electrons. The total spin will be S = 3 if the spins 
are all parallel. But since by Pauli's principle we cannot 
have more than one electron with the same orbital and 
spin quantum numbers to/ and m s and we need to keep 
the angular momenta opposite to the spin, we must have 
a total angular momentum of Ml = m i = — 3. So, 
the maximum orbital angular momentum in this case is 
L = 3. The total J = L - S = in this case. We de- 
note this ground state multiplet term by 2S+1 Lj or in 
this case: 7 F . As usual in the spectroscopic notation, 
one here replaces L = 0, 1, 2, 3 . . . by S, P, D, F . . . . On 
the other hand, a Gd 3+ ion or a Eu 2+ ion would have 7 
/ electrons (half-filled shell) and lead to a 8 S , 7/2 ground 



state. 

Now, when we place the ion in a solid environment, we 
need to worry about how the symmetry breaking of the 
surrounding ions will split these atomic levels. While for 
d electrons, this splitting is strong and leads to so-called 
quenching of the orbital angular momentum, this is not 
the case for 4/ electrons. As mentioned earlier, they are 
shielded from the surrounding ions and stay atomic-like. 
The theory of how the crystal field splits the multiplet 
terms was largely worked out by RacafPH i n terms of 
group theory. 

This gives a very brief summary of the basic atomic 
concepts of / electrons in free atoms and in isolated ions 
as would occur for example for a RE impurity in an ionic 
solid. These atomic multiplet splittings largely deter- 
mine the optical and paramagnetic properties of RE ions. 
While optical transitions between / electrons are strictly 
speaking dipole forbidden in the free atom, the small ad- 
mixture with neighboring ligand orbitals and the result- 
ing symmetry breaking for these ions in a solid environ- 
ment allows weak optical transitions between the ground 
state and the various excited state multiplets. These op- 
tical transitions result in sharp luminescence lines largely 
unaffected by the host and form the basis of many solid- 
state lasers, such as Nd:YAG lasers etc. They are also of 
great interest for RE impurities in semiconductors and 
for optical amplifiers. Because wide band gap semicon- 
ductors like GaN have an advantage for the excitation of 
the RE optical transitions, it has spawned a good deal of 
work on RE-doped GaNpl 

As far as paramagnetic centers, the magnetic dipole 
moment of an ion is determined by the expectation value 
of L + g S in the particular ground state \LSJMj) result- 
ing from Hund's rules. As is well known from textbooks, 
this is determined by the Wigner-Eckardt theorem and 
leads to the Lande g factor. 



B. Early band structure work. 

The earliest approach to RE-V from a band structure 
point of view is to treat the 4/ electrons as core elec- 
trons. In a spin-density functional method, this leads to 
a different potential for spin-up and spin-down valence 
electrons. To the best of our knowledge, the first ap- 
plications o f this approach to RE-V were by Hasegawa 
and YanasePQIZU They successfully predicted already the 
essential elements of the band structure. The valence 
bands are formed from the group V-p orbitals, and the 
conduction bands by the RE-d bands. In all pnictides, 
except for the nitrides, the RE-d band dips below the va- 
lence band maximum (VBM), which itself is located at 
r, at the X-point in the Brillouin zone. These materials 
are thus semimetals, except that the nitrides were found 
to have almost zero gap, so possibly semiconductors, in 
view of the usual underestimation of the band gap by 
the local density approximation (LDA). The drawback 
of this method is that it does not show the 4/ states as 



bands at all. Their location for the occupied spin is in 
fact above the N-2s bands but they are nonetheless not 
included in the band structure picture. The location of 
the empty spin-down states is not calculated. It is as if 
the 4/ electrons live in a different world. 

The same approach was also used by Petukhov et alf^l 
in a rather complete survey of the pnictides, including 
the nitrides. Detailed studies were made of ErAs using 
this approach, in particular because Shubnikov-de-Haas 
measurements were available for the Fermi surface d 74 ! 
The approach was quite successful in determining the 
Fermi surface properties, in particular when spin-orbit 
coupling was included. 

However, it does not accurately describe the spin- 
splitting of the valence bands. If we call the majority 
spin of the / electrons spin-up, then one finds as ex- 
pected that the KE-d spin-up states are also lower than 
the spin-down. This results from the interatomic p — d 
coupling. However, one also finds in these calculations 
that the group V-p VBM also is higher for spin-down 
than spin-up states. This is because the difference be- 
tween the spin-up and spin-down potential simply acts 
more or less uniformly on all electrons in this approach 
and there is no explicit hybridization of the RE-4/ with 
the V-p electrons. In approaches that do take this in- 
teraction into account explicitly, as we will see in later 
sections, the N-p spin-up electrons are pushed up by their 
interaction with the RE spin-up 4/ electrons lying below 
them. The N-p spin-up states thus form antibonding 
combinations with the RE-4/ spin-up states. The spin- 
down states of the N-p, on the other hand, are pushed 
down by the empty 4/ levels above them. Thus, the or- 
dering of the VBM spin states is opposite to that of the 
RE-/ electrons. This has a profound effect on the band 
gaps in the nitrides, as we will see below. 

Let us here mention that the ErAs Fermi surface was 
recently recalculated using a much more involved dynam- 
ical mean field theory,^ which reproduced the Fermi sur- 
face only slightly more accurately, but now also obtains 
the 4/ levels in the band structure picture in the right 
place. 

The main difficulty with 4/ electrons, as already men- 
tioned, is their strongly localized character. As such, 
the local spin density approximation (LSD A), which is 
derived from a free electron gas is no longer valid. If 
one were to do a straightforward LSDA calculation, one 
would find the / levels at the Fermi level, except for 
the Gd case where the half filling will result in a large 
enough splitting to push the levels away from the Fermi 
level even in LSDA. Leaving the core levels out of the 
picture altogether avoids this problem. 
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C. Self-interaction correction (SIC) 



1. Background on the SIC theory 

In Hartree-Fock theory, the terms representing the 
Coulomb interaction of an electron in a given orbital 
with itself are exactly cancelled by the corresponding ex- 
change term. In approximate schemes such as LDA, this 
is no longer the case, and thus there is a residual self- 
interaction error. The self-interaction correction (SIC) 
approach to LDA was first introduced by Perdew and 
Zungerl^Sl It subtracts the self-energy correction for each 
orbital as follows 
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The SIC formalism however has some unusual proper- 
ties. The Hamiltonians determining the one-electron 
states are in principle different for each orbital. So, the 
eigenstates are no longer automatically orthogonal be- 
cause they derive from different Hamiltonians. While 
for a finite atomic system it is straightforward to ap- 
ply the correction, in an extended system with only ex- 
tended orbitals, the correction can be shown to go to zero. 
However, one may find solutions which are localized, i.e., 
break the periodic symmetry, by searching for solutions 
different from Bloch functions with an explicit energy 
minimization approach. This approach was introduced 
by Svane and Gunnarssor l 77 ^ and was first applied to a 
RE system in a study of Ce by Svane.^ 

SIC is somewhat cumbersome because it departs from 
the usual Bloch band structure picture. On the other 
hand, Heaton et alP^ introduced the idea of a unified 
Hamiltonian which, using a projector technique, restored 
a Hamiltonian applicable to periodic Bloch sums of such 
localized solutions on each site. This technique was com- 
bined with the linearized muffin-tin orbital (LMTO) ap- 
proach and further simplified by Temmerman et alpD 



2. Applications of SIC to RENs 

The SIC ap proach was applie d to various RE ele- 
ments: PrJ3l hep GdP and CePHH It was applied 
to Ce chalcogenidesj^ Eu-chalcogenides and pnictidesp^ 
and Yb-compounds (including YbN)P3S3The particular 
choices of Ce, Eu and Yb were because these RE ele- 
ments exhibit a competition between different valencies 
and one of the strengths of SIC is that it can determine 
which f N configuration has the lowest energy. 

The SIC approach was applied to the entire series of 
RENs by Aerts et alPS The paper focuses on determining 
self-consistently how many / electrons need to be treated 
as localized states and hence on determining the effective 
valency. Their results clearly show that all RENs prefer 
the trivalent over the divalent states. As shown in Fig- 
ure [3j there is a systematic trend of the divalent minus 
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FIG. 3. Energy difference between divalent and trivalent RE 
ions in RENs, from Ref. EU 



trivalent total energy: it decreases in the first part of 
the series with a minimum at Eu and then jumps to its 
maximum value in Gd and decreases again toward Yb. 
We will later see that this potential competition between 
trivalent and divalent behavior for Eu and Yb will also 
show up in other approaches. 

This study nicely reproduced the decreasing trend in 
lattice constant with atomic number throughout the se- 
ries, with the exception of Ce. The authors explain the 
trends in bonding and emphasize the hybridization ef- 
fects of the RE-/ with the N-p states. 

In terms of magnetic properties, they find that a small 
magnetic moment opposite to that of the RE is induced 
on the N in the first half of the series and a moment with 
equal sign is induced in the second half. This follows 
the same trend as the RE-4/ orbital moments, which 
are opposite to the spin moment according to Hund's 
third rule. The total spin magnetic moments are integer 
with the exceptions of ErN, TmN and YbN. This implies 
insulating (Tb, Dy, Ho) behavior or half-metallic (Pr to 
Gd) behavior. CeN is found to be non-magnetic and 
metallic. We note that they find the lighter RENs all to 
have a zero gap for majority spin. The minority spin gaps 
gradually increase from PrN (0.53 eV) to EuN (1.46 eV), 
then decrease to 1.11 eV in GdN. The next three, TbN, 
DyN and HoN, are found to have a small gap in both spin 
channels with the minimum gap ranging from 0.05 eV in 
TbN to 0.11 eV in TbN and 0.24 eV in DyN. The last 
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three, ErN, TmN and YbN, are found to have zero gap. 
All of these results are based on densities of states. No 
band structure dispersions are reported in this work. 

As an example, in TbN, Aerts et al. find the occu- 
pied majority spin 4/ bands to lie below the N-2s band, 
while the minority spin filled 4/ bands lie just above it. 
The empty 4/ bands lie just above the conduction band 
minimum (CBM). This differs from the LDA+C/ results 
discussed in Section [lVD which generally find all the oc- 
cupied 4/ bands above the N-2s bands. This indicates a 
very strong effect of the SIC on pushing the occupied 4/ 
bands down. The splitting between occupied and empty 
/ bands of the same spin also appears to be significantly 
larger than in the LDA+C/ results discussed in Sec. |IVD| 

In a closely related paper, Svane et alP^ also studied 
some of the other RE-V and chalcogenides with the SIC 
technique, in particular, the Sm compounds. In the case 
of SmS, SmSe and SmTe, they find a transition from 
trivalent to divalent behavior under pressure. In other 
words, depending on the pressure, the number of / elec- 
trons behaving as localized or delocalized may change. In 
that paper, they also applied the Hubbard-I approaches 
to the spectral functions, which allows one to see the 
atomic multiplet splittings of the 4/ electrons. Unfortu- 
nately, they did not apply this technique to the nitrides. 
We will discuss this approach in more detail in Sec. |IV G| 



D. LDA+t/ 

1. Theoretical background on LDA+U 

An alternative approach for dealing with the localized 
nature of / electrons is the so-called LDA+C/ approach, 
introduced by Anisimov et alP^HsHl l n this approach, the 
orbital dependence of the Coulomb and exchange inter- 
actions for a set of localized orbitals is grafted onto the 
DFT-LDA framework by adding Hubbard-like terms to 
the Hamiltonian and total energy functional: 



Elda+u = Elda + Eif — E, 



dc 



(2) 



The Hubbard terms E\j depend on the occupation num- 
bers of the localized orbitals, or, in a more general for- 
mulation, on their density matrix. In the most general 
open shell case, they are written in terms of Slater's 
Coulomb and exchange integrals^ and Clebsch-Gordan 
coefficients and correspond to the configuration-averaged 
Hartree-Fock expression for the Coulomb energy but us- 
ing an empirically determined screened direct Coulomb 
interation Fq — U . The exchange terms are usually taken 
unscreened and written in terms of one effective exchange 
interaction J. The total energy is then minimized not 
only as a function of the spatial density function n(r) 
but separately as a function of the occupation numbers 
for specific orbitals n ma , or more generally their density 
matrix p a mm i ■ This leads to an orbital dependent poten- 
tial. 



However, because these Coulomb interactions are al- 
ready included in an approximate way in the LDA, 
a double counting correction must be added. There 
has been considerab le discussion of this double counting 
correctio n 1 91 ! 92 ! 94 ^ and there are several slightly differ- 
ent variants of the LDA+C7 approach. One of these is 
the around mean field (AMF) approach while the other 
is the fully localized limit (FLL). The latter is clearly the 
most logical approach for strongly localized states such 
as / electrons. It assumes that in the atomic limit of in- 
teger occupations, the LSDA and LSDA+C/ will give the 
same total energy. 

While in the earliest formulation, the spin effects are 
supposed to result completely from the U and J terms, 
the most common approach is to start already from a 
spin-polarized LSDA and then the U and J terms merely 
need to add the orbital dependence of the Coulomb and 
exchange interactions. The latter are in principle depen- 
dent on the m quantum numbers of the orbitals in an 
open-shell case. This can thus in principle lead to orbital 
as well as spin ordering. Sometimes one spherically aver- 
ages away all these effects,^ while in the most complete 
formulation of Liechtenstein et alP^ these orbital effects 
are included in detail. 

However, the theory still remains a single Slater de- 
terminant Hartree-Fock like theory and thus does not 
deal explicitly with the multiplet splittings discussed in 
the introduction. One may go beyond the Hartree-Fock 
solution for the localized electrons in a so-called dynami- 
cal mean field the ory (D MFT) framework P This will be 
discussed in Sec. IV G| The advantage of the LSDA+C/ 
approach is that it fits more easily into a standard band 
structure approach than the SIC approach. 

In a broad sense, the effect is pretty similar to that 
of the SIC potential, in that it tends to push occupied 
states to lower energy and empty states to higher energy. 
However, unlike SIC where the strength of this interac- 
tion is determined self-consistently within the approach 
in LDA+C/ it is determined by an adjustable parame- 
ter U. In the original work this screened interaction U 
was determined independently from impurity like con- 
strained DFT calculations.^ Recently, a closely related 
linear response approach was proposed by Cococcioni and 
de GironcoliP For the most part, however, it has been 
customary to treat the U parameters as adjustable pa- 
rameters and to either study the behavior as a function of 
U or to determine the values of U based on experimental 
inputs. 

In a formal sense, the relation to a Hartree-Fock treat- 
ment for the localized orbitals with a screened Coulomb 
interaction makes the method closely related to the GW 
theory and also to the recently developed screened ex- 
change and hybrid functionals. 
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2. Applications of LSDA+U to RENs 

Because the application of LSDA+[/ to RENs, other 
RE-V and chalcogenides is straightforward, several 
groups have applied this method. In spite of the sim- 
ilarity of the approach used by these different groups, 
there are significant differences in the results. As usual, 
the devil is in the details. 

One of the earlier attempts to use this method by 
LambrechtP^ did not yet use a full-fledged LSDA+t/ for- 
mulation but rather added energy shifts directly to the 
diagonal Hamiltonian matrix elements in the LMTO ba- 
sis set. The effect on the band structure is the same as in 
the LSDA+C/ method but without the underlying justify- 
ing theoretical framework. In this work, the focus was on 
the question whether ScN and GdN were semiconductors 
or semimetals and shifts were added for both the / and d 
states. These shifts of the / levels were based on the ex- 
perimental X-ray photoemission spectroscopy (XPS) and 
Bremsstrahlung isoch rom at spectroscopy (BIS) for the / 
levels in Gd-pnictidespSS For the d shifts the GW theory 
was used as a guide. It would predict an inverse propor- 
tionality of the shifts to the dielectric constant. The lat- 
ter can itself be obtained from the interband transitions 
in the band structure, so the circle can be closed. As a 
starting point, one needs to know the shift in one mate- 
rial, and ScAs was used for that purpose. That material 
is a semimetal and, as already mentioned above, detailed 
knowledge on the size of the Fermi surface, related to 
the anion-p metal-d band overlap allowed the author to 
determine the required shifts in ScAs. This work success- 
fully concluded that ScN was a semiconductor and was 
the first to provide predictions of the red shift of the gap 
in GdN due to spin polarization. However, this shift ap- 
proach was not sufficiently flexible, too heavily based on 
experimental input and restricted to half-filled / shells. 

We begin our discussion of the LSDA+[/ results with 
the most comprehensive study of the RENs using this 
approach, by Larson et alP^ This paper systematically 
explores the band structures of the entire REN series 
and discusses the magnetic moments. The emphasis of 
the paper is on the question of how the / electron shell 
is occupied. In fact, as the LSDA+J7 formalism for open 
shells can lead to orbital ordering, one must determine 
self-consistcntly not only the number of localized / elec- 
trons (as in the SIC approach), but specifically for which 
m-quantum numbers of the / orbitals the states are oc- 
cupied or empty. In other words, one must make the / 
electron density matrix self-consistent. The trouble here 
is that in principle there can be multiple minima and 
one must find the lowest energy minimum. In particular, 
the authors discuss two separate plausible starting points 
and find which of the two gives lower energy In the first 
one, one assumes cubic symmetry is fully maintained. 
In that case, the energy minimization is dominated by 
the desire to move all the occupied states well below the 
Fermi level. This depends on the specific filling and on 
the crystal field splitting of the / electrons in the octahe- 
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FIG. 4. Top: Energy difference between cubic and Hund's 
rule solutions to the LSDA+17 results. Bottom: Magnetic 
moments in the RENs decomposed into spin and orbital con- 
tribution. From Ref. [29] 



dral environment. In the second approach, one assumes 
a slight modification of Hund's rules. As in Hund's rules, 
the spin is maximized first, but then instead of maximiz- 
ing L, one maximizes L z since the moments are supposed 
to stay fully parallel to each other in a periodic ferromag- 
netic solution. This Hund's rule solution explicitly breaks 
the cubic symmetry by the orbital polarization. It was 
found that the Hund's rule solution had lower energy in 
all cases, except possibly for EuN and YbN as is shown in 
Figure |4| The difficulty with EuN and YbN, as already 
mentioned in the SIC section, is that in these cases the 
divalent solution is close by in energy and during self- 
consistent iterations, the solution may evolve toward an 
erroneous local minimum. 

In contrast to the SIC calculations, the RE-4/ bands 
are found to be always between the N-2s and N-2p bands 
and a gap is obtained in most RENs. Of course, these 
results depend on the choice of U parameters. In fact, 
the origin of the gap in these calculations, as opposed 
to an almost zero gap, results from the inclusion of both 
a Ud and Uf parameter. While the Ut is well justified 
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by the localized character of the / electrons, the Ud re- 
quires some separate discussion. The Uf parameters de- 
termine the splitting of the empty and filled /-states. In 
the Larson paper these were determined rather carefully 
based on experimental dataP2 For GdN, they are based 
on the experimental data on this splitting from X-ray 
photoemission and inverse photoemission for the en tire 
series of Gd pnictides, (GdP, GdAs, GdSb, GdBi).™ 
They also give good results for the occupied / states in 
GdN. The values for the other RENs were then based on 
the assumption that the U / should scale like the atomic 
Coulomb Slater parameter Fq and exhibit similar screen- 
ing. In view of the overall similar electronic structure of 
the outer shells this is an excellent approximation and at 
the same time it incorporates the correct atomic trends. 

Within this approach, the same authors had earlier 
determined that GdN would still be a semimetal with 
t he G d-rf band at X dipping below the N-p like VBM at 
pESI] To account for the experimental observation of a 
gap in GdN, they added a Ud which shifts up the empty 
states relative to the filled states. This parameter was 
adjusted to the best available data at the time for the 
spin averaged gap. It does lead to a prediction of the fer- 
romagnetic red shift of the gap, or the difference in gap 
between majority and minority spins. Although the RE- 
d like bands are not really narrow localized bands, they 
in fact show strong band dispersion. This extension of 
the LDA+J7 formalism to open up gaps should be viewed 
as an empirical correction for the well-known LDA un- 
derestimate of band gaps in most semiconductors. 

As an example, the band structure of GdN is shown 
in Figure [5] The top panel shows the bands occupied 
and empty 4/ bands in red. The lower panels show the 
details along T — X with and without the inclusion of the 
U d - 

Basically, the importance of Ref. [55] is that it clearly 
establishes that Hund's rules stay valid even for periodic 
RE compounds. This is born out in practice by the fact 
that the RE in RE-V definitely have an orbital moment 
and hence break the cubic symmetry. The slight twist 
on Hund's rules of making it compatible with a periodic 
ferromagnetic solution, however, results in an interesting 
situation of approximately zero magnetic moment for Sm 
rather than for Eu. As we will see in Section [V] this is 
confirmed by the experiments for SmN in the ferromag- 
netic state. 

Among the RENs, GdN and the related pnictides have 
received more attention than any of the others. In part, 
this is because these are the easiest to deal with theoret- 
ically because they have a half-filled shell so even LSDA 
gives reasonable results. The LSDA+/7 ap proach w as ap- 
plied to t hese compounds by Dua n et al.llfl SHlQSl Ghosh 
et al.pa Larson and Lambrechtp 9 ' 101 ' 1 " 6 ' Mit ra and 
LambrechpMoI] and Abdelouahed and AlouaniPMini 
We will discuss the differences between these papers in 
later sections focused on the magnetism (Section [v]) and 
the optical properties (Section VI). First, we complete 
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FIG. 5. Band structure of GdN in LSDA+J7. Top panel, from 
Ref. 1291 The /-bands are colored red, the blue and purple in- 
dicate mixed / and other orbitals character. Solid and dashed 
lines indicate majority and minority spin. Second and third 
panel, from Ref. 11011 Detail of the bands near the gap with 
and without Ud- 



with the correlation of the / electrons. 



E. Hartree-Fock, Hybrid functionals 

As was emphasized in the two previous sections, the 
main error of LSDA for / electrons is the self-interaction 
error. By treating the / electron Coulomb terms sepa- 
rately at the Hartree-Fock level, LSDA+C/ corrects the 
self-interaction error, and SIC of course explicitly re- 
moves it. The drawback of LSDA+fJ is that it is some- 
what empirical and ad-hoc as the / Coulomb interac- 
tions are just added onto the theory. One might want to 
use a more rigorous first-principles approach which treats 
all electron orbitals on the same footing. Hartree-Fock 
theory is self-interaction free but misses correlation en- 
tirely. It is known to strongly overestimate band gaps 
in semiconductors, and to tend to overestimate spin- 
splittings and magnetic moments. Its application to / 
electron systems has been delayed by technical difficul- 
ties. Nonetheless, recently, it has become possible to ap- 
ply Hartree-Fock to periodic solids with Gaussian basis 
sets and r ecen tly also to include / electrons in the CRYS- 
TAL code.Eil 

This also opened the way toward hybrid functionals. 
The hybrid functional approach, i.e., mixing some frac- 
tion of Hartree-Fock with LDA, has long been popular 
in the quantum chemistry community but has recently 
gained considerable ground in solid state physics ap- 
plications. Two functionals have been commonly used, 
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B3LYFEMXn] and more re cently the HSE 115 ' 116 ' (Heyd- 
Scuseria-Ernzerhof) functional. The latter mixes in 
about 25 % exact exhange, but cuts off the long range 
part of the exact exchange, and thereby uses effectively a 
screened exchange. This screening is done with an error- 
function type cut-off. The method was recently imple- 
mented i n the V ASP (Vienna Ab-initio Simulation Pack- 
age) cod c! 117 * 1 -^ and has since gained a lot of popularity 
as it was found to be promising to significantly improve 
the band gaps. It has not yet been widely applied to / 
electron systems. B3LYP and Ha rtree-Fock were both 
applied to GdN by Doll ct alP^ using the CRYSTAL 
code. Very recently Schlipf et alP^fll implemented the 
HSE functional within an FLAPW code and applied it 
to GdN. 

As expected, Hartree-Fock gives a large band gap of 
about 5 eV in GdN and also overesti mat es the spin- 
splitting of the / electrons. Doll et alP-^ find the oc- 
cupied / levels at about 15 eV below the VBM and the 
empty spin-down levels at about 22 eV above the VBM. 
Interestingly, the B3LYP gave different solutions with 
somewhat different total energy. The band structures in 
the lowest two of these solutions show a strange extra 
band near the Fermi level with rather different disper- 
sion to the well-known N-p or Gd-<i like bands. In one 
case, it occurs in the majority spin and in the other in the 
minority spin. In a third solution, which however, they 
find to have a 2.4 eV higher total energy, they find an in- 
sulating band structure with the usual indirect (T — X) 
gap of about 0.7 eV for majority spin. The empty and 
occupied / levels lie about ±7 eV relative to E F and are 
in reasonable agreement with experiment and LSDA+J7 
calculations. The origin of the strange additional solu- 
tions in B3LYP is not clear. 

We should mention one other Hartree-Fock plu s cor- 
relation energy study of GdN by Kalvoda et alPS3 This 
paper however focuses entirely on the cohesive energy and 
does not give any information on the band structures. It 
treated the / elect rons as core states. 

Schlipf et alP^I used the HSE functional and obtained 
results very close to the LSDA+[/ results. They find the 
/-levels at about ±6 eV relative to the VBM. They find 
a very nearly zero band gap, with a slight positive gap 
of 0.01 eV at the experimental lattice constant (4.988 A) 
and a slight negative gap of —0.06 eV at their theoretical 
equilibrium lattice constant (4.963 A). Thus, like Duan 
et al. li£SJ they obtain a solution at the brink of a metal- 
insulator transition. Their smallest direct gap at X is 
0.90 eV for the ferromagnetic state and 1.17 eV (average 
of spin-up and spin-down gap) above Tq- This red-shift 
of th e ga p of 0.27 eV is quite close to the results of Larson 
et alPB or Trodahl et al. 64 of about 0.4 eV. 



F. The GW method 

The most accurate and rigorously first-principles com- 
putational method available today for band structures 



is Hedin's GW approach!^!! Here G and W stand for 
the one-electron Green's function and screened Coulomb 
interaction respecively. While the band structure eigen- 
values of the Kohn-Sham equations in density functional 
theory are in principle only Lagrange parameters related 
to the one-electron wavefunctions' orthogonality and nor- 
malization, used in minimizing the total energy as a func- 
tion of density, they do not strictly speaking represent 
excitation energies. The theory for quasiparticle excita- 
tions, i.e., extracting or adding an electron to the system, 
is conceptually the correct theory for addressing photoc- 
mission and inverse photoemission theory. Within this 
many-body theoretical framework, GW is the first order 
approximation for the self-energy operator in an expan- 
sion of the screened Coulomb interaction. The equations 
of the GW method can be summarized as follows. The 
quasiparticle excitation energies Ei are given by 



1 



+ / d s r'^ xc {v,T' ,Ei)^i{r') = Eii,i{r), 



(3) 



where vsr(r) is the Hartree potential representing the 
classical interaction with the nucleus and all electrons, 
S xc (r,r',w) is the self-energy operator which is non- 
local and energy dependent and ipi ( r ) is the quasiparticle 
wavefunction. The self-energy operator is given by 

£ xc (r,r» = < kSG (T,i>,u-u')W(T,i' ) w')e- iSu ' 

(4) 

in which Go is the one-electron Green's function of 
the corresponding independent particle approximation, 
which is usually the Kohn-Sham equation, 
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V 2 +v H (r) +v xc (r) 



^(r)=e^(r), (5) 



and with the screened Coulomb interaction, 

W(t,t',oj) = J e- 1 (r 1 r'\oj)v{r" 1 r)d 3 r" (6) 



in which s is the dynamic dielectric response function and 
v the bare Coulomb interaction (l/|r — r'|). Schemati- 
cally, one can write W = (1 — uHo) _1 w in which Ho is the 
independent particle polarizability, which itself is given 
in terms of the Green's functions as 



no(l,2) = -iGo(l,2)Go(2,l). 



(7) 



Here we used a short-hand notation where 1 stands for 
{ri,<7i,ii}, i.e. position, spin and time of particle 1. In 
practice, these equations are Fourier transformed over 
time and lead to a convolution in the frequency domain 
and, instead of directly using the position dependent 
functions, the corresponding operators are expanded in 
suitable basis functions. 
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While this formalism has been available since Hedin's 
work of 1965, applications to real systems have lagged 
because of the considerable technical challenges in im- 
plementing it, in particular because the frequency depen- 
dent dielectric function is needed. The first applications 
in the 1980s were restricted by pseudopotentials and ap- 
plied only to standard semiconductors but led to much 
more accurate band gaps than LDA. Implementing it to 
systems with localize d d and / electrons proved a fur- 
ther challenge II22H125I Typically, the equations are solved 
only in first order perturbation theory starting from the 
LDA, which one calls a "one shot Go Wo" approximation. 
In the most recent incarnation, which is called the quasi- 
particle self-consistent GW (or QSGW) approach, a new 
exchange correlation potential is extracted from the T, xc 
and leads to new one-electron Green's functions, which in 
turn lead to a new Yi xc . These equations are then solved 
self-consistently and the quasiparticle energies Ei then 
become equal to the £j. 

This method was implemented in the FP-LMTO 
framework by van Schilfgaarde et al ! 126 ' 127 ' It is impor- 
tant to note that this is an all-electron method with- 
out the need for pseudopotential approximations, so that 
core-valence exchange is properly included. It has al- 
ready been shown to give accurate band gaps for all 
standard semiconductors and some transition metal com- 
pounds. The method was applied to / el ectro n systems 
including GdN and ErAs by Chantis et alPU 

The QSGW method gives a band structure that sup- 
ports the results of prior LSDA+C7 calculations in general 
terms. That is, the filled and empty / orbitals are moved 
far away from the band gap region. Looking in more 
detail at the gap, the direct application of the QSGW 
method gives a gap of 0.22 eV for GdN. Because in most 
semiconductors, QSGW is found to overestimate the gaps 
slightly, a hybrid approach between LDA and QSGW us- 
ing 80% QSGW and 20 % LDA as exchange correlation 
was used and gives a gap of only 0.05 eV. This method 
thus once again gives a gap very close to the metal insu- 
lator transition in GdN. The direct gaps at X are 0.46 eV 
for majority and 1.48 eV for minority spin and thus give 
an average of 0.97 eV, close to the experimental absorp- 
tion edgeP Rec ent m easurements place this gap slightly 
higher at 1.3 e V 55 * 64 * which suggests that the uncorrected 
QSGW may be closer to the experiment in this case than 
with the correction factor of 0.8. So clearly GW does 
quite well on the band gaps and also supports the HSE 
calculations discussed in the previous section. 

The main problem with the QSGW approach at 
present is that it appears to overestimate the majority 
and minority spin splitting of the / levels. In partic- 
ular, the empty / levels lie significantly higher than in 
the LSDA+f/ and HSE results. This is also obvious for 
other 4/ syst ems like ErAs and metallic Gd as shown by 
Chantis et alP^H These authors blame this discrepancy, 
as well as the remaining overestimate of the QSGW gaps 
in most cases, to the underestimate of the screening by 
the random phase approximation used in Eq. [7| On the 



other hand, this seems to not affect the quality of the 
energy levels near the gap. 



G. Dynamical Mean Field Theory 

In all previous methods, the / electrons are still taken 
into account at the Hartree-Fock level. In other words, 
in a many-body theory framework, they correspond to a 
single determinant solution. These methods thus do not 
do justice to the intricate splittings of the / levels one ob- 
serves in atomic systems and also in solid state systems 
for any but the simplest half-filled case. The problem 
is that a multi-electron framework in which the multi- 
plet splittings appear seems conceptually totally sepa- 
rate from a one-electron band structure picture. The 
dynamical mean field theory (DMFT) provides a way 
to overcome this difficulty. The key is to realize that 
what a band structure really represents is the one-particle 
excitation spectrum. One may then calculate the one- 
particle excitations between different multi-electron con- 
figurations (or multiplet terms, to use the terminology 
established in the introductory section) corresponding to 
specific quantum numbers L,S,Ml,M$ for the isolated 
/ electrons first. These define a well defined one-particle 
Green's function for the isolated / systems without in- 
teraction with the other electrons in the system. They 
define an on-site self-energy for the / electrons. The 
DMFT approximation then consists in assuming that 
the self-energy is not k-dependent and can be solved by 
a so-called impurity solverP^H The interaction with the 
remaining electrons can then be switched on assuming 
these local interactions occur at each site and this way the 
broadening effects on the / levels due to their interaction 
with other (band electrons) in the system are included as 
well as the indirect effect of the / electrons on the rest 
of the system. Typically, the quantities compared with 
experiment in this case are the spectral functions which 
can be directly compared to photoemission and inverse 
photoemission. 

The DMFT has many applications beyond 4/ elec- 
tron systems, including any system where more complex 
many-body effects occur in the localized electron system. 
The dynamical aspects of the problem are important for 
example in Kondo systems, where they can give rise to 
both upper and lower Hubbard bands, but also a peak 
at the Fermi level. It is als o im portant for the actinides, 
as shown by Savrasov et alP^I For 4/ electrons its main 
importance is that it provides a way to put the multiplet 
splittings and band structures on the same footing in a 
single comprehensive picture. This particular application 
of DMFT was previously introduced by Lichtenstein et 
alpSl and called the Hubbard-I approximation. It goes a 
step beyond LSDA+C/ but is in the same general spirit of 
including the relevant atomic physics, but grafted onto a 
DFT formalism. 

The DMFT approa ch w as applied to Ce-pnictides by 
Laegsgaard and SvaneP^ The approach was recently ap- 
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plied to ErAs as already mentioned and led to results 
in good agreement with the Fermi surface Shubnikov- 
de Haas experiments.^ At the same time it gives quite 
different /-level splittings than LSDA+J7. Experimen- 
tally, from X-ray absorption and emission spectra for 
edges, i.e., 3c? — > 4/ transitions, it is clear that the 
spectral shapes are indeed dominated by these atomic 
multiplet effects. 

Interesting hybridization effects of the localized multi- 
plet split levels with other bands were observed in p ho- 
toemission spectra of EuNi 2 P2 by Danzenbacher et alP^I 
These effects are also typical for heavy fermion systems 
involving light RE such as Ce. In the LSDA+t/ band 
structure studies by Larson et al.p^EuN stood out as a 
possibly interesting case from this perspective. In fact, it 
was found that this case was one of the exceptions to the 
Hund's rules in which a solution of the equations main- 
taining the cubic symmetry had lower energy. This how- 
ever, turns out to be due to the proximity for Eu of the 
divalent and trivalent solutions. This proble m w as pre- 
viously pointed out by Johannes and Picket j n t, ne 
work of Larson et al.pS it led to a prediction of a rather 
unusual band stucture with a band of mixed d and / 
character which changes its character as a function of k 
in the Brillouin zone crossing the Fermi level. EuN was 
thus targeted as one of the RENs worthy of a more de- 
tailed investigation. Recently this came to fruition in a 
paper on EuN combining various theoretical approaches 
compared to various X-ray spectroscopies.^ It was found 
that the LSDA+C/ fails for this system, while GW gives 
a correct prediction of the existence of a band gap but 
only the DMFT in the Hubbard-I approximation can sat- 
isfactorily describe the details of the / electrons and their 
hybridization effects on the Eu d bands. 

Even though the DMFT comes closest to combining 
the atomic multiplet aspects of the / electrons with the 
band structure picture, it is still limited in that it is 
geared toward studying the one-electron excitation spec- 
tra, as occurs most directly in photoemission and in- 
verse photoemission, and to some extent in X-ray ab- 
sorption and emission. In other words, it can treat exci- 
tations in which the number of /-electrons changes by 1: 
jN jN-i Qr jN _^ yiv-K However, it does not yet 
allow a proper treatment of the optical excitations. In 
that case, one is often interested in excitations between 
multiplet terms within a single f N configuration. These 
types of excitations are mostly important for RE impuri- 
ties but may also become important to fully understand 
the optical properties of RENs. 

Recently, a constrained DFT approach was proposed 
by Hourahine^ 3 -^ which could potentially deal with such 
situations. It minimizes the DFT total energy under the 
constraint of a particular set of L, S, Ml, Ms quantum 
numbers. In other words, it constrains the energy op- 
timization by constraining the expectation values of L 2 , 
5 2 , L z and S z . It has not yet been applied to RENs and 
was only implemented in a tight-binding type of band 
structure approach. For the sake of completeness, how- 



ever, we mention it here as an approach which has poten- 
tial to provide further progress on this difficult problem. 
It is in some sens e a generalization of earlier work by 
Eriksson et alP^ to incorporate orbital polarization so 
as to impose Hund's rules, also discussed by Solovyev et 
alPSl 

Having done a tour of the main approaches to treat 
the electronic structure of the 4/ electrons within a band 
structure framework and their applications to RENs, we 
will discuss in the following sections some details on how 
these methods give us results to be compared to measure- 
ments concerning magnetism and interband transitions. 
But first we complete this section with a discussion of 
lattice dynamics. 

H. Pressure effects and phase transitions 

The sensitivity of the GdN band structure to the lattice 
constants was already pointed out by var ious authors and 
in particular emphasized by Duan et alP^l As we men- 
tioned earlier, GW and HSE also find GdN to be at the 
brink of a metal insulator transition. Hence, studies as 
a function of pressure or for films under tensile stress, as 
can be realized in epitaxial layers, are important. Ab- 
delouahed and AlouanPSH took this a step further and 
investigated the possibility of phase transitions between 
different crystal structures. They show that under tensile 
stress there would be a transition to zinc blende GdN. It 
is somewhat doubtful that such a transition can ever be 
reached because it requires a negative pressure, but it 
is nonetheless of interest as a hypothetical compound in 
the context of Gd-dop ed G aN. On the other hand, Ab- 
delouahed and AlouanP^II predict also a phase transition 
to the wurtzite structure under applied pressure. This is 
rather surprising given that the wurtzite phase has lower 
coordination that the rocksalt phase and remains to be 
confirmed by experiment or other calculations. They pre- 
dict a transition pressure of 68.3 GPa in GGA+J7 and 19 
GPA in GGA. 



V. MAGNETIC STUDIES 

The REN series was originally a goldmine for neutron 
scientists, who found ferromagnetic materials with a sim- 
ple cubic structure very easy to investigate through neu- 
tron scattering. This outcome was in contrast to work on 
the pure RE metals which have more complex crystallo- 
graphic structures. Through such studies, the ferromag- 
netic nature of most of the RENs had been found already 
in the early 1960s!^2l The easy axis of magnetisation, as 
controlled by the intrinsic anisotropy of the RE, was also 
determined. In a few cases, the lack of single crystals left 
some uncertainties in the magnetic structures, because 
powder samples only give the direction of the ordere d mo- 
ments with respect to the scattering vector. ^ 7 * 138 ! How- 
ever, the general findings were in agreement with theo- 
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retical predictions based on an atomic approach to the 
RE magnetic properties, resulting from crystal field cal- 
culations performed in the ground state multiplet J 

Since less than ten years ago the magnetic properties 
of the RENs have been reinvestigated owing to the com- 
bined progress in the growth of thin films and in ab initio 
calculation techniques. The RENs are indeed now con- 
sidered as model systems for band calculation and used 
to check the ability of different models to describe the lo- 



calised RE states (see Section IV I . The central question 



of the exchange mechanism in the RENs has also been 
addressed by these theoretical calculations. 

In the following we will give some examples of the more 
recent studies on the magnetic properties of REN thin 
films, which can be considered as real advances after so 
many years of limited progress. First, in Section |V A| 
we recall the experimental tools which are used to es- 
tablish magnetic properties, with emphasis on the tech- 
niques for thin films. Next, we give details about some 
of the more significant data, obtained predominantly on 
epitaxial thin films and mostly on GdN (Section VB|. 
A great step forward in understanding the magnetism 
of the RENs was taken when the prominent role of Vjv 
as carrier dopants was experimentally demonstrated. In 
Section |"V C| we describe the recent experimental results 
obtained on other RENs, in particular SmN and EuN. Fi- 
nally, in Section |"VD| we discuss the different theoretical 
models that have been proposed to explain the magnetic 
properties of the RENs. 



A. Experimental details 

Standard magnetometry with a Superconducting 
Quantum Interference Device (SQUID) or Vibrating 
Sample Magnetometer (VSM) is the first step needed 
to determine the magnetic properties of the REN films, 
and these techniques have been the basis of all exper- 
imental studies. The magnetic moment, the magnetic 
ordering temperature and the magnetic anisotropy can 
all be derived, however the diamagnetic or paramag- 
netic contributions from the substrate or the buffer and 
cap layers may introduce some uncertainties inherent 
in thin films. For weakly anisotropic fcrromagnets like 
GdN, the correction for parasitic magnetic signals is easy, 
while for the other RENs with non-collinear ferromag- 
netic arrangements or strong anisotropics, the correction 
can be problematic. 15 Having recourse to complemen- 
tary techniques that avoid the substrate contribution en- 
tirely is becoming the rule, mainly Ferromagnetic Res- 
onance (FMR) and X-ray Magnetic Circular Dichroism 
(XMCD). In addition, each of these techniques provides 
its own specific information. 

FMR brings valuable information about both the mag- 
netisation and the anisotropy constants, even on poly- 
crystalline films, through the application of a static mag- 
netic field along different directions with respect to the 
film plane! 14 " ' 141 ' 



X-ray absorption spectroscopy (XAS) is becoming an 
essential technique for probing the magnetic properties 
of materials, in particular with XMCD. The element se- 
lectivity of XMCD ensures that only the magnetic ele- 
ment of interest will be analysed. The energy selectivity 
permits access to the magnetic polarisation of each elec- 
tronic level separately; 4/ and 5d for the RE and 2p for 
the N. There are two energy ranges, which probe differ- 
ent aspects of the RE magnetism. The M^s absorption 
edges in the soft X-ray spectral range probe transitions 
to the 4/ orbital empty states. Therefore, the XMCD 
difference signal intensity is proportional to the 4/ mag- 
netic moment. With the help of specific sum rules, it 
is in principle possible to separate the spin an d orbital 
contributions to the magnetic moment.' 142 * 143 ' The 1/2,3 
absorption edges in the hard X-ray spectral range probe 
the transitions to the 5d band empty states, which may 
be polarised through the intra-atomic f-d exchange inter- 
action. In metallic RE compounds the 5d states are con- 
duction electron states and contribute to the exchange 
mechanism between the localised 4/ spins. In the case of 
the RENs, the L-edge XMCD signal provides indirect in- 
formation about the 4/ magnetic moment and also about 
the other sources of polarisation in the lattice. In addi- 
tion, the RE XAS peak positions are well separated for 
different valence states, for both spectral edges. The en- 
ergy separation is important in the case of anomalous 
REs which may exhibit mixed valence states (Ce, Sm, 
Eu, Tm, Yb). We will show in Section VC an illustra- 
tion of this situation in the case of EuN. 

The existence of magnetic polarisation at the N site 
can be obtained independently by measuring the XMCD 
signal at the N K-edge. The spectra involve the Is — > 2p 
electronic transition and probe the 2p character of the un- 
occupied states. Leuenberger et al. observed a very large 
XMCD signal at the N site in GdN The corresponding 
hysteresis loop from SQUID measurements exactly re- 
flects the Gd 4/ magnetisation derived from the Gd M 4 
edge XMCD signal. The same authors have investigated 
the XMCD spectra at the Gd £2,3 edges in the case of a 
strained film of GdN (see details in Section VB3)p3l 



A full XMCD investigation of EuN at the M 4j5 and 
-^2,3 edges of the Eu ion has given a fundamental set of 
informatio n 48 ' 49 ' which we will discuss later. Theoretical 
calculations have been of great support to understand the 
experimental results and have in turn benefitted from the 
data obtained by measurements. 

Antonov et al. have provided an in-depth theoret ical 
study of XMCD in GdN for various absorption edgesP^ 
They used an ASA-LMTO Hamiltonian and applied both 
LDA and LDA+C/ methods. In addition, they investi- 
gated the effects of the core-hole on the spectra and on 
the surface states. While they find good general agree- 
ment in the main features for all edges, they point out 
the need to include multiplet splitting effects to fully un- 
derstand the high energy shoulders in the fine structure 
of the M 4i5 edges. The L 2 ,3 edges are explained in detail 
in terms of the various spin-orbit split states. In addi- 
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tion, Antonov et al. study the effect of electric quadrupole 
and magnetic dipole transitions on these spectra. Ab- 
delouahed a nd A louani have also studied the GdN i 2 ,3 
edge spectraPH For the N K-edge, the XMCD clearly 
requires the inclusion of core-hole effects and even then 
Antonov et al. could not fully explain the fine st ruct ure 
observed in the experiments of Leuenberger et alP^H Fi- 
nally, the surface states make a sizable contribution to 
the first peak at 400 eV in the N K-edge of GdN. 
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FIG. 6. The measured coercive field against the inverse of 
crystallite volume in GdN thin films. From Ref. [53] 



In the following we will give an overview of the main 
parameters that show a significant impact on the mag- 
netic properties of the RENs. We show how a better 
control of these parameters in the recent studies on thin 
films has permitted some advance in the understanding 
of the mechanism of exchange interactions in the RENs. 
These data all focus on GdN as the model material. 
The compound GdN was indeed the first to be reinvesti- 
gated by several groups. There soon appeared some dif- 
ferences in the results published mainl y by two groups, 
in GottingerP^ and in Wellington! 53 1 65 1 More recently, a 
group in Kobe has reported on the magnetic properties of 
GdN thin films 55 and another group in Cambridge has re- 
ported peculiar properties due to the existence of another 
Gd-N phase depending on the preparation conditions.^ 
Other important differences between these studies con- 
cern the transport properties, which will be described in 
the dedicated Section IVT1 below. 



1. Control of the impurity 

In the older preparations of bulk material, a level of 
oxygen contamination was unavoidable. Even in thin 
films where the growth atmosphere can be controlled very 
accurately in ultra-high vacuum systems, a small amount 
of contamination cannot be excluded. It is well known 
that O levels below 1% cannot be detected through ion 
beam analysis techniques such as Rutherford backseat - 
tering or secondary ion mass spectroscopies. However, 
thanks to these older experimental reports, the effect of 
oxygen contamination on the magnetism has been anal- 
ysed and understood. The ferromagnetic ordering is can- 
celled with 5% O in the latt ice and short range antifer- 
romagnetic order sets inP^ 6 -! Therefore, any increase of 
Curie temperature Tc cannot be accounted for by the 
presence of oxygen impurities. 

2. Crystallite size 

GdN is a soft ferromagnet, with the low temperature 
in-plane coercive fields of recent thin films reported as 
ranging from approximately 10 Oe to 220 Oe! 47 | 53 | 55 | 59 | (i5 1 



Senapati et al. have reported an enhanced coercive field 
of several hundred Oe in GdN grown by reactive dc mag- 
netron sputtering at high power, as well as the presence 
of a small exchange bia s in sa mples that may containg a 
strained phase of GdNp3D32] For films with the relaxed 
GdN structure, Figure [6] shows a linear dependence of 
the coercive field on crystallite size, suggesting that the 
magnetic reversal mechanism may take place through nu- 
cleation of domains at defect sites and that the intrinsic 
magnetic anisotropy of defect-free GdN is near 15 Oe.^ 

3. Strain effects 

In 2005, Duan et a l 11221 investigated the effect of strain 
on the band structure of GdN with a LSDA+ U calcu- 
lation process. They demonstrated that the system ex- 
hibits a half-metallic band structure at the equilibrium 
lattice constant and that a semimetallic or semiconduct- 
ing character develops with increasing lattice constant. 
They show that the magnetic properties are extremely 
sensitive to the volume variation. The calculated magni- 
tude of the exchange parameter is reduced with increas- 
ing volume, and so is Tc- 

There are a few experimental works which report on 
the effect of an increased volume on the magnetic proper- 
ties of GdN. Leuenberger et al. investigated the influe nce 
of a lattice expansion on the Gd £2,3 XMCD spectraP^Sl 
Tc is significantly reduced (30 K instead of 60 K), sug- 
gesting a reduction of the exchange interaction. The film 
deposited by N + plasma-assisted reactive ion beam sput- 
tering at room temperature shows a 9% larger unit cell 
volume with respect to the bulk parameter when the film 
is deposited at 450° C. It is shown that the 5d states are 
polarised and carry a magnetic moment. Apart from 
this reduction in Tc, the authors observe another inter- 
esting effect. The XMCD signal amplitude ratio IL3/L2I 
= 3 is unexpectedly high in the volume expanded film, 
while statistically a ratio of 1 is expected and observed 
in most Gd systems, being imposed by the degeneracy of 
the 2pi/ 2 and 2p 3 / 2 core states. The authors attribute 
this effect to a reduction of the L2 signal amplitude. The 
influence of the lattice constant enlargement is not called 
on directly, but a more fundamental reason is proposed. 
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P(N 2 ), Torr 

FIG. 7. Curie temperature Tc and conductivity a just above 
Tc vs N2 growth pressure. The error bars result from the 
less than 100% contribution from the Gd ion moments to the 
moment in the paramagnetic phase. From Ref. 1131 



A change in the 5d polarisation is due to a change in 
the relative occupation of the spin up and spin down 5d 
bands, predominantly of majority spin Gd 5d ti g char- 
acter. The L2 edge mainly probes the Gd 5d%/2 and the 
5d tig states have a stronger j—3/2 character. XMCD 
is thus a very accurate technique to give reliable infor- 
mation about the polarisation of the 5d Gd states, and 
consequently about the effect of changes in the exchange 
interactions. As a comparison, Preston et al. measured 
the 1/2,3 XAS spectra of a GdN (001) epitaxial fil m wi th 
a bulk-like lattice parameter and a Tc of 70 KP^H A 
1^3/^2! ratio of 1 was measured as expected, which pro- 
vides another proof that the 9% volume expansion of the 
Leuenberger et al. sample affected the exchange interac- 
tion. Preston et al. claim that the anomalies detected in 
the L 2 edge XAS and XMCD amplitudes are related to 
a modification of the 5d ti g states, which are involved in 
the bottom of the conduction band and experience ex- 
change splitting below Tc ■ This is in agreement with the 
theoretical calculation by Duan et alP^S 

The Wellington grou p al so investigated a volume ex- 
panded sample of GdNpSD The sample was a polycrys- 
talline film grown by ion beam assisted deposition which 
exhibited a lattice parameter of 5.12 A, corresponding to 
a 2.4% tensile strain. Tc was drastically reduced down 
to 20 K and the magnetisation measured in-plane satu- 
rated only above 4 T, suggesting some disorder-induced 
anisotropy. The grain size was indeed quite reduced 
(3 nm) compared to the microstructure of the bulk- 
like films (10 nm) reported earlier.^ This was confirmed 
thr ough FMR experiments, performed on both types of 
£i m LL40] j n k th films, a strong uniaxial anisotropy per- 
pendicular to the film plane was found, which is at- 
tributed to strain effects. The magnetoelastic anisotropy 
contribution was smaller in the lattice expanded sample. 
The origin of this contribution is not yet understood. 

In summary, the strong strain effect observed in the 



magnetic properties of GdN shows that the exchange in- 
teraction is sensitive to the relative occupation of the 
spin states in the 5c? band, and that the polarisation of 
the conduction band is indeed involved in the Gd 4/-4/ 
exchange, suggesting a carrier-mediated exchange inter- 
action brought about by doping. We will discuss this 
issue in the next section. 

4- Stoichiometry - Effect of N vacancies 

A general agreement about the Curie temperature of 
GdN has been found at near 70 K (see Table [Tl]) , but all 
authors of such studies agree that the RE/N stoichiome- 
try plays a major role. Similar to the effect of unwanted 
oxygen cont aminat ion, a lack of nitrogen was shown to 
decrease j^ EHEini Very recently, Plank et alP^ carried 
out a systematic study of N deficient polycrystalline thin 
films of GdN, resulting in the important discovery that 
a N deficit enhances the average Tc* values and is equiv- 
alent to carrier doping. These carriers may be polarised 
by the localised magnetic moments and participate in an 
exchange interaction process of the RKKY type. 

In the earliest study of the Wellington group it was 
reported that the rate of N 2 flow during growth and the 
sample conductivity were intimately linked, suggesting 
that the importance of the right stoichiometry is the key 
to understanding the transport and magnetic properties 
of GdNP*l It was only in the last couple of years that 
a systematic study of the Vn became possible, owing to 
progress in optimising the preparation conditions and a 
strict control of the nitrogen partial pressure.^ The re- 
sults of Figure [7] delivered through the influence of the 
growth pressure are spectacular. It was found that an 
increase of the conductivity by eight orders of magni- 
tude occurred by reducing the N 2 growth pressure from 
2.6xl0 -4 Torr to 7xl0~ 6 Torr, resulting also in a near 
tripling of the Tc value to ^200 K. This also demon- 
strates that some hope exists to utilise GdN as an ef- 
fective material for spin-polarising devices by driving Tc 
even higher through proper management of the Vjv con- 
centration. 

More experimental work has established a link between 
the magn etic be haviour and the Vn concentration in the 
same way! 14 * 14 ^ The authors claim the presence of an- 
other phase which may be antiferromagnetic, but this 
conclusion needs some more experimental work. 

C. Magnetic properties of other RENs 

Whilst the magnetic behaviour of the reference com- 
pound GdN may be near a full understanding, some of 
the other RENs have not even been reinvestigated since 
the 1970s. In recent years significant progress has been 
made towards understanding a few heavy RENs (DyN 
and ErN) and light RENs, in particular two special mem- 
bers of the series, SmN and EuN. The results are all de- 
rived from thin films, either polycrystalline or epitaxial. 
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TABLE II. Structural and magnetic properties of RENs. FM = ferromagnetic, AFM = antiferromagnetic, PM = paramagnetic, 
VV PM = Van Vleck paramagnetic. Moment is the maximum value of the moment at low temperatures and /x e ff is the 
paramagnetic moment derived from a Curie- Weiss law fit above 8 P . 
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a Hulliger 1979 ^ and refs. therein 
b Hulliger 1978 Sand refs. therein 
c Vogt and Mattenb erger 1993^1 and refs. therein 
d Preston et al. 2007 G^H 
c Meyer et al. 2008 EH 
f Ruck et al. 2011 1331 
s Xiao and Chien 1996 ESU 
h Li et al. 199TE21 
1 Osgood III et al. 1998 E3] 
j Plank et al. 2011 EU 
k Nakagawa et al. 200^23 

I Leuenberger etal. 2005^ 
m Si et al. 200'£2>H 

II Scarpulla et al. 2009 
° Natali et al. 2010 EU 

P Senapati et al. 2011 ED 
1 Thiede et al. 2011 
r Wachter et al. 1998 EH1 
s Yamamoto et al. 20o4±H 
* Nakagawa et al. 200d23 
u Meyer et al. 2010 ESS 
v Degiorgi et al. 199(f2] 
w Kasuya and Li 199$^ 



In agreement with previous studies performed on bulk 
polycrystalsp^l semiconducting DyN in thin film form 
was demons trat ed to undergo a ferromagnetic transition 
below 25 KpS In the case of ErN thin filmsPH the T c 
of 6 K is also in agreement with previous investigations, 



however more details have been established owing to the 
(111) texture of the films. The magnetisation curves of 
the textured thin films measured in an applied magnetic 
field parallel and perpendicular to the plane of the films 
are not compatible with the original assumption that 




T(K) 

FIG. 8. Temperature dependence of the magnetisation of 
SmN after cooling without (ZFC) and with (FC) an applied 
field of 5 kOe. The inset shows the result of a fit of the FC 
susceptibility (x = M/ H) in a Van Vleck approach. From 
Ref . [T5l 



the Er magnetic moments would be perpendicular to the 
(111) direction. Neither is the magnetisation compatible 
with the moments being aligned along the (111) direc- 
tion, which would be in agre ement with calculations of 
the single ion anisotropyP^I Instead, the magnetisation 
mostly agrees with an alignment of the magnetic mo- 
ments along the (100) direction, similar to the magn etic 
moment arrangements found in HoN and DyNPBSS 

Among the light RENs (RE=Ce, Pr, Nd, Sm, Eu), 
SmN and EuN are of special interest because of the prox- 
imity of excited states. SmN involves the Sm 3+ ion, 
which in the free ion state ( 6 H 5 / 2 ) presents a magnetic 
moment fi = gJ = 0.71 /i^. However, due to the small 
energy separation between the ground and lowest excited 
state multiplets, the angular momentum J is in most 
cases not a good quantum number and the magnetic mo- 
ment has to be defined as /i = — (< L z > +2 < S z >) fiB- 
The spin and orbital angular momenta are equal in mag- 
nitude with opposite sign, so that they cancel each other 
giving rise to a zero total magnetic moment. A zero to- 
tal moment is expected in SmN, and any magnetic order 
is indeed difficult to evidence. Previous experiments on 
bulk SmN were rather in favour of antiferromagnetic or- 
dering instead of a ferrom agnetic on e with a small mag- 
netic moment (see Table |n|) ) 3 I 7 I 15 3 The recent experi- 
ments on thin films definitely show that the compound 
undergoes a ferromagnetic transition below 30 K, shown 
in Figure |8p^ Above this temperature, the magnetisation 
curves and low field susceptibility are in agreement with 
a Van Vleck paramagnetism. In the ordered state, the 
magnetic hysteresis loop exhibits an increasingly large 
coercivity when lowering the temperature while the mag- 
netic moment remains very small. This is reminiscent 
of the well-known behaviour of fcrrimagnetic materials 
close to their compensation point. Due to the polycrys- 
talline structure of the films, there is no saturation of 
the moment at high field because of the huge anisotropy 



20 




I ■ 1 ■ 1 ■ 1 ■ 1 ■ ' — 

20 40 60 80 100 



T(K) 

FIG. 9. Eu Ms-edge amplitudes at 1131 eV (black filled cir- 
cles) and 1134 eV (open red circles) are measures of the mo- 
ments on the Eu 2+ and Eu 3+ ions in EuN, respectively. The 
black line is a fit to the Eu 2+ moment by the Brillouin func- 
tion appropriate for a half-filled 4/ shell in a field of 50 kOe. 
From Ref. \M 
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FIG. 10. Eu L-edge amplitudes vs temperature in EuN. From 
Ref. EU 



of the Sm ion. However, the analysis of the data within 
a ferromagnetic scenario is undisputable. In addition, 
the nearly zero moment is in agreement with theoretical 
calculations!^ The compound SmN is thus of great in- 
terest as a ferromagnetic material without a fringe field, 
as has already been p ointed out for the metallic alloy 
(Sm^Gd^AlaPSE^However, unlike (Sm 1 _ 2; Gd 2: )Al 2 , 
SmN is a semiconductor. 

EuN also involves an irregular RE becauseEu in its 
free ion trivalent state has no magnetic moment (J — 0). 
EuN is thus expected to be a Van Vleck paramagnet, 
with the Zceman interaction mixing several multiplets of 
Eu 3+ owing to the weak separation energy of 460 K be- 
tween the J = ground state and the J = 1 excited state. 
Therefore, a constant magnetic susceptibility is expected 
at low temperature. Old data were never conclusive be- 
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cause of the presence of O impuritiesP Recent experi- 
ments were achieved on epitaxial thin films.^ Unexpect- 
edly, the magnetisation data show a Curie-like tempera- 
ture behaviour of the susceptibility. The explanation is 
provided by XMCD experiments, using a combination of 
the data from two absorption edges of Eu. The M4 5 ab- 
sorption edge directly probes the magnetic state of the 
4/ orbitals, while the L2.3 absorption edge probes the 
empty states of the 5d band, which is expected to be 
polarised through the intra-atomic / — d interaction. A 
small amount of Eu 2+ is detected in the M-edge XAS 
spectra, corresponding to about 1% of the total Eu sig- 
nal. The XMCD spectra are dominated by the divalent 
impurity signal, owing to the intrinsic high paramagnetic 
moment of Eu 2+ (7 fig). Figure [9] shows the temperature 
dependence of the M-edge XMCD signal for the two va- 
lence states, resembling a Brillouin function for Eu 2+ and 
quasi-constant for Eu 3+ . The origin of these two contri- 
butions is found in the L-edge experiments on the same 



sample. Figure 10 reproduces the temperature variation 
of the L2 XMCD spectra. In agreement with the M-edge 
results, the d-orbital divalent contribution strongly de- 
creases with increasing temperature, while the trivalent 
contribution remains practically constant. An ab-initio 
calculation of the XAS and XMCD spectra was under- 
taken to understand this behaviour. The comparison be- 
tween the experimental spectra and the simulated one 
allowed the relationship between the Eu 3+ and the Eu 2+ 
bd signal amplitude to be evaluated. Clear evidence for 
a linear dependence is established, suggesting that part 
of the magnetic polarisation of the 5d orbitals in EuN 
is due to the Eu 2+ impurities through interatomic 5d-5d 
exchange. 

It is interesting to note a small number of early studies 
that were reported on N-containing RE compounds with 
elevated Tc values owing to the inclusion of a significant 
fraction of dopants replacing the N ions. GdNi-^C^ was 
reported to be ferromagnetic with a Tq value of 190 K, 
and material with an additional defi ciency of N was re- 
ported with a T c as high as 340 K . |160 | 1( * 1 1 In addition, 
replacing up to 30% of O in EuO with N augments Tc 
from 69.5 K to as much as 77 KP^I This latter result 
was attributed to a decreasing RE-RE distance with an 
increasing N concentration, resulting in a stronger ferro- 
magnetic exchange interaction between the RE cations. 
These reportsfurther support the possibility of enhanc- 
ing Tc values in the RENs through appropriate control 
over the strain or carrier concentration by wa y of doping, 
similar also to the case of Gd-doping in EuOP^ 3 - 1 



this compound has received the bulk of the recent atten- 
tion. The magnetism of GdN is usually studied together 
with the rest of the Gd pnictide series. 

Before we get into a detailed discussion of the 
LSDA+f/ results, we want to discuss the question of why 
the ferromagnetism of GdN is of special interest. First, 
it is somewhat surprising that the compounds which are 
definitely semimetallic in the Gd pnictide series - GdP, 
GdAs, GdSb, GdBi - are all antiferromagnetic, while 
GdN, the only one which is now believed to be a semi- 
conductor, is ferromagnetic. Usually, one finds antiferro- 
magnetism in insulators and ferromagneti sm i n metals. 
Second, as pointed out by Kasuya and Lijl^ the simi- 
larity of the ferromagnetic properties of GdN to those of 
EuO is surprising. Although both systems have a half 
filled / shell (because Eu in EuO is divalent), the posi- 
tion of the energy levels is quite different. In EuO the 4/ 
bands lie above the 0-2p bands and thus form the VBM, 
while in GdN, as we have seen, the majority spin 4/ 
bands lie well below the N-2p bands. In a simple model 
of the ferromagnetism the exchange interaction between 
nearest neighbour Gd in second order perturbation the- 
ory is given by 



t 2 

t 2 - or df 



(8) 



Here Idf is an intra-atomic d—f exchange coupling, while 
tdf is a hopping parameter from d to / orbitals on neigh- 
bouring sites, S is the spin and AE^f is the splitting of 
the d and / levels. This model thus would predict a much 
smaller exchange coupling and hence smaller Tc in GdN 
than in EuO. It is therefore quite surprising that both 
materials have similar Tc, ~69 K for EuO and 58-70 K 
for GdN. This observation led Kas uya and Li to propose 
a novel exchange coupling scheme^S arising from 4th- 
order perturbation theory and described by 



(4) 



-4/, 



<-(/' 



SAE pf AE 2 



(9) 



D. Theories of Magnetism in RENs 

In this section, we discuss further in detail studies of 
the ferromagnetism in RENs, mostly on GdN because 



involving p — d as well as / — d hopping. Instead of using 
such intricate perturbation coupling schemes, spin den- 
sity functional theory attempts to obtain the exchange 
couplings directly from the total energy calculations, as 
we now describe. 
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TABLE III. Comparison of exchange interactions in Gd-pnictides (GdX) from different papers. 
GdX Ji (meV) J 2 (meV) 





Larsorj^] 


DuarQ 


Mitra 


Lf] 


Larson" 


Dua^ 


Mitra c 


Li d 


GdN 


1.84 


0.86 


0.42 


1.74 


-0.87 


-0.14 


-0.36 




GdP 


1.13 


-0.17 


0.34 


0.60 


-1.35 


-0.74 


-0.82 


-0.92 


GdAs 


0.98 


-0.22 


0.12 


0.22 


-1.52 


-0.91 


-1.03 


-0.95 


GdSb 


0.80 


-0.51 


0.15 


0.14 


-1.95 


-1.13 


-1.22 


-1.63 


GdBi 


0.69 


-0.66 


0.01 


-0.11 


-2.45 


-1.37 


-1.44 


-1.71 



a Larson and LambrechllinUFP-LMTO,!// + U d 
b Duan et al.ES3FLAPW,t/f only 
c Mitra and LambrechlPZlFP-LMTP, U f + U d 
d Li et al.EH expt. 



The common approach used in several of these papers 
is to map a generalised Heisenberg Hamiltonian 

H = -J2 J v S i- S h (10) 

ij 

to the total energies of a series of different magnetic con- 
figurations, for example ferromagnetic (FM) order, anti- 
ferromagnetic ordering along [001] (alternating layers of 
up spin and down spin along a [001] direction, or so-called 
AFM-I), antiferromagnetic ordering along [111] (AFM- 
II), and other hypothetical ordering schemes. Within 
a model with a certain number of parameters, say first 
and second-nearest neighbours, one can then easily write 
down expressions for the total energy in each magnetic 
configuration. Calculating these energy differences from 
first principles using LSDA+J7, one then extracts the 
exchange parameters from the model. From these ex- 
change parameters, the Curie or Neel temperatures are 
estimated using either the mean field approximation or 
some more sophisticated statistical method, for example 
the Monte Carlo approach. 

To demonstrate the convergence of the model, one typ- 
ically tries to show that other configurations beyond the 
ones used to extract the model parame ters are also well 
reproduced. For example, Larson et alP^ consider also 
the AFM-III ordering, which consists of alternating 2 lay- 
ers of up spin with 2 layers of down spin along [001] 
and show that its energy difference from the FM case is 
well reproduced if they determine the nearest and second- 
nearest neighbour parameters J\ and Ji from the AFM- 
I-FM and AFM-II-FM energy differences. Duan et alPH 
include a third-nearest neighbour interaction J 3 instead 
but find it to be small. There are slight variations to 
this approach, as some authors apply the method to a 
quantum Heisenberg Hamiltonian with S = 7/2, while 
others use a classical Heisenberg Hamiltonian. At the 
mean field level these map exactly to each other, so it 
is straightforward to translate the exchange interactions 
from one study to another even if the original papers use 
slightly different definitions. In the discussion below we 
will refer to classical Heisenberg spins and assume the 



sum is over all sites i and j in Eq. |10| i.e., each pair is 
counted twice. 

Before we discuss the differences between the results 
of the different studies, let us point out that, apart from 
GdN, all agree on the lowest ground state for the Gd- 
pnictides to be AFM-II ordering, in agreement with ex- 
periment, and that all find GdN to be ferromagnetic, also 
in agreement with experiment. 

Duan et alP^I emphasized the volume dependence of 
the properties of GdN. They noted that GdN in their ap- 
proach comes out very close to a metal-insulator transi- 
tion and hence carefully studied the volume dependence. 
At the calculated equilibrium lattice constant, they ob- 
tained a semimetallic band structure, but at slightly ex- 
panded lattice constant a band gap was obtained. They 
noted the strong dependence of the exchange interactions 
on the lattice constant. As mentioned above, this sugges- 
tion was picked up by Leuenberger et al.p^who stud- 
ied strained GdN experimentally and indeed found a be- 
haviour close to a metal-insulator transition. 

In a second paper, Duan et alP^H studied the exchange 
interactions up to third nearest neighbour interactions 
in the Gd-pnictide series and explained the increasingly 
stronger antiferromagnetic second nearest neighbor cou- 
pling parameter in terms of superexchange via the in- 
tervening N. They suggested that the nearest neighbour 
exchange interaction is RKKY in nature and obtained a 
ferromagnetic interaction for GdN and antiferromagnetic 
coupling for the other pnictides. 

Larson and Lambrecht^^ independently obtained sim- 
ilar results but used an additional Ud as already men- 
tioned, which turns GdN into a semiconductor even 
without an increase in lattice constant. Mitra and 
LambrechtP^ improved on the results of Larson with a 
more careful method of extracting the exchange inter- 
actions. The key here was to recalculate the ferromag- 
netic/antiferromagnetic energy differences in exactly the 
same unit cell so that systematic computational errors 
cancel. They obtained much better agre ement for the 
exchange interactions than Duan et alP^ A comparison 
between these three calculations is given in Table |III| 
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TABLE IV. Comparison of Neel and Curie- Weiss temperatures obtained by different groups. 
GdX Tjv (or T c in GdN) (K) T C w (K) 



Lars 


o 


43 


Duar 


^ Mitra [] 


if] 


Larson" 


Duan 6 


Mitra c 


Li d 


GdN 


47 


c 




34 [ 


8 e 


58 


67 


37 


11 


81 


GdP 


22 


g 




12 


13 9 


15.9 


21 


-25 


-3' 


4.0 


GdAs 


24 




14 


17 


18.7 


10 


-31 


-18 


-11.8 


GdSb 


31 




20 


20 


23.4 


-8 


-50 


-22 


-31.3 


GdBi 


39 




22 


23 


25.8 


-25 


-62 


-32 


-34.0 



a Larson and LambrechlfinUFP-LMTO,^/ + U d 

b Duan et al.ESlFLAPW JJf only 

c Mitra and LambrechlPZlFP-LMTP, U f + U d 

d Li et al.E3 expt. 

c 0-7T cw 

f Monte Carlo simulation 



Although different codes were used by the Lam- 
brecht and Duan groups, respectively a full-potential lin- 
earized muffin-tin orbital (FP-LMTO) method and a full- 
potential linearized augmented plane wave (FLAPW) 
method, this is not the main origin of the differences. 
The band structure methods give essentially converged 
results for a given choice of Hamiltonian. The remain- 
ing differences for GdN in particular can thus be clearly 
traced to the use of a Ud parameter opening up the band 
gap, as used by the Lambrecht group. 

Besides the Neel temperature T/v, one can also calcu- 
late the Curie- Weiss temperature Tew within the mean 
field model and compare both to the experiments. Ar- 
guably, the calculations of Mitra and Lambrecht obtain 
the trend in both T/v and Tew slightly more accurately. 
The predictions of the different groups for the critical 



temperatures and Tew are compared in Table IV It is 



important to note here that the experimental values for 
the exhange parameters were in fact obtained by Li et 
a i 03] f rom experimental values for Tew arid Tjv analysed 
within mean field theory. The exchange parameters are 
found using 



Tw = -J-2/k Bl 
T CW = (4Ji + 2J 2 )/k B 



(11) 



respectively within the mean field approximation and 
thus determine directly J 2 and J%. As the Tew values 
are obtained from the high temperature susceptibility, 
the mean field theory is better justified for them than 
for T/v. The results of Duan et alLLMI f rom both Monte 
Carlo and mean field methods show that the actual crit- 
ical temperatures are typically underestimated by about 
30 % by the mean field theory. Overall, one can say 
that the magnetic exchange interactions in the series of 
Gd pnictides is fairly well understood on the basis of 
the LSDA+L 7 method. Similarly, the same approach was 



also used for the Eu-chalcogenides^Sl and in that system 
also reproduced the trend in T/v and Tew in good agree- 
ment with experiment. Again, one finds a switch from 
ferromagnetic to antiferromagnetic AFM-II behaviour in 
the series, although the theory predicts the switch to oc- 
cur between EuO and EuS, while in the experiment the 
switch occurs between EuS and EuSe. 

Ghosh et alP^ calculated magnetic moments of the Gd 
monopnictides but not the exchange interactions. The 
magnetic moments also show small differences between 
t he d ifferent groups. For example, for GdN, Ghosh et 
al.[i05| report a spin (orbital) momen t of 6.79 (0.13) [i B 
on Gd and -0.14 [i B on N. Duan et alPS report 7.038 fj, B 
as the total magnetic moment on Gd and Larson et alP^ 
give an / contribution to the Gd spin moment of 6.93 /is, 
a d contribution to the spin moment of 0.081 [ i B a nd a N 
moment of —0.083 fi B . Mitra and Lambrechti^ report 
only Gd-d and N induced magnetic moments of 0.08 and 
—0.08 hb, respectively. Of course, all these calculations 
result in approximately 7 [i B per Gd. However, different 
authors have emphasised different ways of decomposing 
the total magnetic moment. In principle, one expects a 
zero orbital moment for Gd 3+ in. agreement with Larson 
et alpsD Mitra and LambrechtP^ emphasized the exact 
cancellation of the induced magnetic moments of the Gd- 
d orbitals with the N moment in the ferromagnetic but 
not the antiferromagnetic state. They thus coined the 
description of GdN as an "antiferromagnet in disguise" . 
In fact, the Gd-d moments are surrounded by N moments 
in perfect antiferromagnetic alignment, and vice- versa. 

In fact, the basic picture that emerges from all these 
calculations is that the exchange interactions in GdN 
arise from the small induced moments on Gd-d and N- 
p orbitals which are coupled to the large 4/ moment on 
Gd through intra-atomic exchange coupling. The 4/ mo- 
ments themselves are too localised to have any direct cou- 
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pling between each other. There are different ways to de- 
scribe the effective exchange couplings depending on how 
one maps these on an effective Hamiltonian. The com- 
mon approach is to map them on a classical Heisenberg 
model with spins only on the Gd atoms. As already men- 
tioned, it is within this approach that Li and Kasuya^^J 
introduced a phenomenological model involving various 
orders of perturbation theory and which is also used in 
the LSDA+C/ work discussed above. Unfortunately, no 
direct correlation between this perturbation theory and 
the LSDA+C/ calculations has been found, in particu- 
lar not for the high order th eory proposed for GdN. On 
the other hand, Duan et alf^Sl find some perturbation 
theory justification for the trend with anion diameter in 
J2, which they describe in terms of a superexchange via 
the anions. The larger the anion, the larger the overlap 
and hence the increasing trend in magnitude of the Ji 
with increasing anion size, which is in fact observed in 
all LSDA calculations. Unfortunately, the trends in J\ 
are less well understood and there is poorer agreement 
between the different calculations of J\. 

However, one may also calculate the exchange interac- 
tions in a model including spins on both Gd and N. In 
fact, Mitra and LambrechtP^ analyse this different ap- 
proach using the exchange interactions calc ulated using 
Liechtenstein's linear response approach Because this 
model is evaluated within the atomic sphere approxima- 
tion (ASA) to LMTO, it includes empty spheres on the 
interstitial sites (introduced to obtain a better space fill- 
ing with spherical potentials). In many ways, this is a 
more direct approach. It calculates the exchange inter- 
actions directly from the Green's functions for a reference 
magnetic configuration, usually the ferromagnetic ground 
state. It considers small transverse fluctuations from this 
model to define the exchange interactions, rather than 
the total energy differences between different collinear 
spin arrangements. This model makes no a priori as- 
sumptions on what is the correct Hamiltonian or which 
spins to include. The analysis of Mitra and LambrechlP^ 
shows that this model justifies the simpler approach in- 
cluding only Gd spins for all Gd pnicitides other than the 
nitride. For the nitride, one finds exchange interactions 
of similar magnitude between Gd, N and empty spheres. 
Thus, they suggest, one perhaps should view GdN as a 
ferrimagnet rather than a ferromagnet. 

Interestingly, these different viewpoints lead to simi- 
lar conclusions about the critical temperature in GdN, 
namely that it is very small. Within a mean field ap- 
proach to the ferrimagnetic model (not incuded in that 
paper but worked out later by the author), one obtains 
19 K as the critical temparature, whereas the spin on 
Gd only model gave 11 K in Mitra and Lambrecht's 
work. While the improved computational details even- 
tually gave a convergence between the different group's 
results for the exchange interactions in the other Gd pnic- 
tides, the differences became more obvious for GdN and 
in fact worsened the agreement with experiment. The 
estimate of LarsorP^ gave a Tc of 67 K, in good agree- 



ment with the experiments which gave values typically 
around 70 K, although Li et alPI place T c at 58 K. 
U sing a more accurate Monte Carlo approach, Duan et 
alPlobtained T c =34 K. Even if we apply a Monte Carlo 
vs. mean field reduction factor of 0.7 to Larson's results 
we obtain Tc=47 K. In Mitra's results, we obtain only 
11 K in mean field or 8 K if we include the Monte Carlo 
reduction factor. Clearly the difference between the cal- 
culations of Mitra and Duan must result from the fact 
that in Duan's case GdN has a semimetallic band struc- 
ture while in Mitra's case the band structure is semi- 
conducting. However, we know that a semiconducting 
state is in agreement with experiments! 64 ! 65 ! While an 
RKKY carrier-mediated explanation for the ferromag- 
netic nearest neighbour interaction makes some sense in 
a semimetal, it does not in a strictly insulating material. 
Even so, it is clear that the interactions must come from 
the d electrons. One can thus think of a sort of virtual 
excitation to the d like conduction band which then me- 
diates the interaction between the localised 4/ moments. 
This is pretty much the visi on i ncorporated in Kasuya 
and Li's perturbation theoryP^Il \y e are faced with the 
problem that improving the methodology led to better 
agreement with experiments for the other pnictides but 
worse agreeme nt f or the nitride! For this reason Mitra 
and LambrechtP^ also checked the case of Gd metal with 
the linear response method and arrived at Tc=285 K, 
(including the 0.7 correction factor) in quite good agree- 
ment with the experimental value of Tc— 297 K. This led 
the authors to the conclusion that pure GdN must have 
a much lower Tc and that the disagreements with exper- 
iment might result from extrinsic factors such as Vjv. 

A first atte mpt at studying the latter was reported 
by Punya et alP^S This study showed that Vm indeed 
increases the nearest neighbour interaction and reduces 
the second-nearest neighbour antiferromagnetic interac- 
tion, leading overall to an increase in Tc by about a factor 
2. This, however, does not yet take into account the pos- 
sibility of longer range, more itinerant interactions from 
the presence of free carriers in the Gd d-band. 

Sharma and Nolting^^ 1 present a separate study of the 
free carrier contribution to the exchange interactions in 
GdN. Their many-body Green's function formalism is 
somewhat intricate and its relation to the Liechtenstein 
linear response formalism is not entirely clear. For ex- 
ample, simply shifting the Fermi level in the latter ap- 
proach to simulate free carriers did not lead to a strong 
increase in TcrP^ In Sharma's approach the Gd / elec- 
trons are represented as localised spins only and the band 
structure of the remaining electrons is extracted from an 
LMTO Hamiltonian without explicit inclusion of the / 
electrons P^SJ The coupling to the localised spins is then 
treated with a more advanced Green's function method 
than the simple linear response. Nonetheless, this work 
presents some evidence that the discrepancy between the- 
ory for perfect insulating GdN and experiment may in- 
deed stem from the carrier-mediated interactions. 

The result of these studies is to bring us almost full 
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circle to the asserti ons from the early work on GdN by 
Wachter and Kaldis, 1 ^ who claimed that "semiconduct- 
ing GdN (n/Gd<C 1CP 3 ), if it were possible, would be 
an antiferromagnet in zero field. However, the lowest 
carrier concentration obtained until now makes GdN a 
semimetal" . While the recent theory work on pure GdN 
which best agrees with the known experimental data, 
namely that it is in fact a semiconductor, still gives a 
ferromagnetic rather than antiferromagnetic state, it is 
one with a very low Tc- 

Evidence for carrier mediated ferromagnetism and 
even a second antiferromagnetic phase was recently re- 
ported by Senapati et alP^ in films with Tc=60 K. On 
the other hand, in another recent experimental studyj^a 
much lower Tc of 37 K was reported for GdN. This value 
was based on using Arrott plots and on fairly thin films 
of 30-95 nm thick. They attributed the differences to the 
other recent studies by Granville et alP^and Leuenberger 
et alP^ to differences in stoichiometry, film thickness and 
grain size. 

In a recent hybrid functional calculation (using the 
HSE functiona l 115 * 116 |), somewhat higher Tc values of 
55 K in the mean field approximation an d 42 K in a ran- 
dom phase approximation were reportedP^I So- perhaps 
there is finally some convergence between theory and ex- 
periment on the value of Tc- The exchange interactions 
Ji and J2 as defined earlier in the HSE calculation are 
found to be 1.09 eV and 0.17 eV respectively. So, in 
contrast to the LSDA+t/ results, they find a positive J2 
and somewhat larger J\, but we do not yet know how 
effective HSE is for the other Gd pnictides or EuO. 

Other magnetic properties are also of interest and in 
this context it is important to mention the magnetocrys- 
talline anisotropy energy (MAE) determined in the recent 
study of Abdelouahed and AlouanfSS] f or Gd, GdN and 
GdFe2- They find a different easy axis in fee Gd (along 
[111]) than in GdN (along [001]). They find the right or- 
der of magnitude for the MAE in hep Gd and show that 
GGA+U, GGA and GGA with open core calculations 
lead to results that differ by orders of magnitude. Al- 
though they did not extract explicitly the anisotropy pa- 
rameters, the difference between the hard and easy axes 
appears to be of order 30 /ieV in hep Gd, 5 /ieV in fee 
Gd and only 0.5 ^eV in GdN. They thus predict a rather 
large tunability of the magnetic anisotropy in GdN s with 
various amounts of N. 

As for the other RENs, a systematic study of ferromag- 
netic versus antiferromagnetic ordering has not yet been 
done. Johannes and PicketflSS studied the exchange cou- 
plings for EuN and EuP. Their perspective is that Eu 3+ 
is a nominally non-magnetic ion with J = 0, as discussed 
earlier. Their LSDA+[/ band structure is similar to the 
ones recently obtained in the GW approximation!^ In 
agreement with the SIC calculations of Horne et al.p^ 
they find a trivalent solution, but in disagreement with 
Larson et al.p^ they find a way to satisfy Hund's rules 
with maximum L z . As in the case for the other RENs, 
they find a semimetallic band structure for EuN when no 



explicit corrections are made to the RE-d bands. They 
find the ferromagnetic energy to lie 4.9 meV below the 
AFM-II configuration. This is comparable to the results 
for GdN ob tain ed by LarsonPSH The paper by Johannes 
and Picketl^ 3 - 3 - on EuN and EuP includes an interest- 
ing discussion of the magnetic exchange interactions, in- 
cluding a more complex effective Hamiltonian than the 
Heisenberg Hamiltonian, which explicitly includes spin- 
orbit coupling terms. 



VI. ELECTRONIC AND OPTICAL 
PROPERTIES 

It is in the determination of the electronic behaviour, 
the electronic band structures, where the recent thin- 
film studies of the RENs have made the most important 
advances. As discussed in the introduction the reason- 
able agreement about the magnetic states of most RENs 
found in 20th century studies was in direct contrast to 
the uncertainties about their electronic band structures 
and transport properties. In Section |IV| the most recent 
band structure calculations were presented, and here we 
review the experimental evidence for those band struc- 
tures and for the conducting states; are they semiconduc- 
tors or semimetals? In view of both ambient-temperature 
electronic and low-temperature spintronic potential, it is 
of interest to answer these questions for both the mag- 
netically ordered and the paramagnetic phases, and to 
establish the effects of magnetic order on the transport 
properties. Thus in this section we review more recent 
transport studies probing the states near the Fermi level 
and optical and X-ray spectroscopy investigations used to 
probe the band gaps and states deeper into the valence 
and conduction bands. 

Historically there was an enormous variation in trans- 
port measurements. It was realised that defects played a 
defining role, which itself is the signature of a semicon- 
ductor. Furthermore, essentially all results found in early 
studies^sH showed low-energy, presumably free-carrier 
absorption and an enhanced absorption at higher energy 
interpreted as signalling interband edges in the 0.7-2.9 eV 
range. There was substantial disagreement even among 
measurements on any one of the RENs as regards the 
band gaps as well as the density of free carriers. The 
gaps reported are shown in Table [Vj based on an assort- 
ment of experimental works. 

In 1990 there was a thorough investigation of YbNpl 
suggested at the time to promise a heavy fermion 
state. The results were interpreted to indicate YbN is 
a semimetal with an optical gap of 0.1 eV, now with a 
low carrier concentration absorbing relatively weakly at 
lower energies. The study in this case was extended to 
well above the gap, reporting several higher-energy inter- 
band transitions. Even to date this is the only of these 
materials to have anything more than an optical band gap 
reported. It is interesting that the same group reported 
metallic conduction in TbN, with an absorption onset 
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TABLE V. Interband edges in the paramagnetic phase 
{T>T C ) 
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a Hulliger 1979 ^ and refs. therein 
b Xiao and Takai 199^201 
c Richter et al. 201& 
Shalaan and Schmidt 
e Trodahl et al. 200^ 
f Yoshitomi et al. 201lEIl 
s Bommelli et al. 199^23 
h Wachter et al. 199^11 
1 Bommelli et al. 199$^ 
j Wachter et al. 199^^ 
k Sclar 1964SU 
1 Azccm et al. 20lJE! 
m Brown et al. 20l£22] 
n DeGiorgi et al. 199(13! 



near 0.75 eV that they interpret as a plasma edgeP^ 



We start in the next section with a review of transport 
studies, which investigate the excitations at low energy 
where the differentiation between semiconductors and 
semimetals is most apparent. Following that discussion 
we review the few optical studies that have been reported, 
mostly to determine optical band gaps. XAS/XES inves- 
tigations of the full band structure, delineating the par- 
tial densities of filled and empty states, are then reviewed 
in Section [VI B | Although these do not have the resolu- 
tion to signal a clear gap between valence and conduction 
bands, they show considerable detail in the overall den- 
sity of states and offer a more complete comparison with 
gross behaviour of the calculated band structures. 
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FIG. 11. The conductance of GdN as a function of the pres- 
sure of N2 during deposition, recorded during a growth under 
a continuously diminishing pressure. From Ref. 1651 



A. Transport and the optical band gap 

1. GdN 

Recent conductivity measurements of GdN still show 
considerable variation among different growth condi- 
tions, but nonetheless the results are within a somewhat 
smaller range than historical data. Reported ambient- 
temperature values for GdN range fr om 10 ficm in 
high-temperature grown epitaxial film J 16 ' 47 ' 53 ' to 1 »cm 
in polycrystalline films grown at lower temperatures li^lfiH 
The dominant dopants are Vn, as has been clearly 
demonstrated by the dependence of the conductivity on 



the partial pressure of N2 during growth (Figure 111. 
It is expected that each vacancy will bind either one 
or two electrons, thus acting as either single or double 
donors at all temperatures, with the carrier concentra- 
tion rising a s the bo und electrons are released at higher 
temperature! 13 * 59 -^ The few Hall measurements that are 
reported confirm that th e carriers are electrons in the 
conduction bandP^HZES] The reaction at the surface of 
deposited Gd is slow, so that the ratio of N 2 /Gd flux 
must exceed 100 to achieve conductivities as low as 10 
Scm" 1 . Such low conductivities suggest carrier concen- 
tration of order 10 18 -10 19 cm^P^It is tempting to assign 
that also to the order of Vat density, though the possible 
existence of traps would prevent accumulation of all do- 
nated electrons in the conduction band. Epitaxial films 
are grown at higher temperature, for which the low for- 
mation energy of Vn, along with the likely reduced con- 
centration of electron traps, ensure they are more heavily 
dope d th an ambient-temperature grown polycrystalline 
films Hall effect data on the more conductive films 
show electron concentrations of order 10 20 -10 21 cm~ 3 , 
with at most a very weak temperature dependence that 
is all but masked by the field dependence of the extraor- 
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dinary Hall effect PES There appear to be no Hall effect 
data on less conductive films, but evidence for carrier 
freeze-out comes from strongly activated conductivities 
in the paramagnetic stated 

The temperature coefficient of resistance (TCR) varies 
from weakly positive to strongly negative as the ambient- 
temperature conductivity falls, signalling thermally ac- 
tivated transport. However, conventional semi-metallic 
behaviour has also been reported in films with high car- 
rier density.^ 3 * 47 * 53 ! A well-established resistive anomaly 
forms a peak at Tq , as seen in Figure |12[ with the resis- 
tivity falling in the ferromagnetic phase by up to an order 
of magnitude P^Ml Both metalli c and semicond uctor signs 
of TCR are reported below r c ! 13 ' 16 ' 47 ' 53 '^Ml The data, 
and especially their sensitivity to the growth temperature 
and N 2 pressure, suggest that GdN is a semiconductor in 
both the paramagnetic and ferromagnetic phases. 

A sizeable magn etores istance has been observed in 
GdN by two groups! 53 * 59 ^ At 100 K, the non-saturating 
magnetoresistance reaches 5% in a field of 50 kOe. As 
shown in Figure |12| with an applied field the anomalous 
peak in the resistivity seen at Tc is reduced and shifted to 
higher temperatures. The result is an enhanced magne- 
toresistance near Tq which approaches 35% in epitaxial 
filmsP 

A clear ambient-temperature (paramagnetic) optical 
absorption ed ge of 1.3 eV in GdN has been reported 
by two groupg^ES in the past few years. Both stud- 
ies also found a red shift of the gap to 0.9-1.0 eV 
deep in the ferromagnetic phase, (Figure 13 1, a result 
of the narrowed majority-spin gap in the ordered spin 
state. Very recently, evidence was found also for a red 
shift of both majority-spin and minority-spin gaps in an 
AIN/GdN/AIN heterostructureP^ It is likely this nar- 
rowing leads to an enhancement of carrier concentration 
as evidenced by the resistive anomaly at Tq. The para- 
magnetic and ferromagnetic gaps have combined to tune 
the value of Ud in a LSDA+[/ band-structure calculation. 
The refined spin-ordered band structure then shows a fi- 
nite but relatively small absolute gap of about 0.4 eV 
between the majority-spin VBM at T and the CBM at 
x¥^ Sub-gap absorption in both optical studies seems 
to be small, implying carrier densites of less than 10 21 
cm~ 3 , in agreement with a doped semiconductor or very 
low-carrier-concentration semimetal. Interestingly, the 
gap is close to those measured in early studies, but in 
enormous disagreement with the 4.1 eV claimed more re- 
cently for GdN on films that had not been protected by 



a passivating cap 
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2. Other RENs 

There appear to be only four other RENs for which 
recent optical or transport measurements have been re- 
ported. An optical gap of 0.95 eV was measured for 
EuNj^l which has been of theoretical interest for its pos- 
sible divalency and for a suggested hidden spin order in 
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FIG. 12. Resistance (normalised to the value at 298 K) of 
a 1000 A GdN film as a function of temperature in differ- 
ent magnetic fields (left panel). The resistivity at 298 K 
is ~10 mfl cm. Right panel: Magnetoresistance (R(H) — 
R(0))/R(0) vs magnetic field at temperatures below, near, 
and above Tc ■ From Ref . [59] 



> 
03 



50 

-50 
-100 



w -150 
-200 



-250 
-300 



-rrpr 



i i i | i i i i | i i i 



_L 



-L 







50 100 150 200 250 300 
Temperature (K) 

FIG. 13. Temperature dependence of the direct band-gap 
energy of 95-nm thick GdN. AE g corresponds to the change 
in the band-gap energy versus the room-temperature value. 
From Ref. 



the more likely trivalent stateP^The existence of the in- 
terband edge and the relatively weak free carrier absorp- 
tion below the gap are strong evidence that EuN is non- 
metallic, in direct disagreement with some of the calcu- 
lated band structures.^ CeN is reported to be metallic,^ 
but in contrast it has also been reported to be a low- 
conductivity semiconductor with an optical gap of 1.76 
e y ll2PJ BoN has been reported to show optical gaps in the 
1.48 eV to 3.32 eV energy range, but with substantial ab- 
sorption below the gap suggesting a very h eavy doping, 
and possibly even metallic carrier densities Most re- 
cently, DyN has been shown to have an optical gap of 
1.2 eV with a blue-shift where Vjy are present associated 
with the Moss-Burstein effect P^H 
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Recent temper atur e-de pend ent resi stivi ties have been 
reported for SinN, 1 ^ DyN^ and ErNP^All three show 
resistivities typical of heavily doped semiconductors, and 
negative TCRs with anomalies at their Tq that are very 
much weaker than found in GdN. 

Sumarising the transport/optical evidence, then, there 
is strong evidence that SmN, GdN, DyN and ErN are 
semiconductors. EuN shows a clear interband edge at 0.9 
eV, with sufficiently weak subgap absorption to suggest 
that it is semiconductor. HoN shows a clear interband 
absorption edge charateristic of a semicondutor but the 
free carrier absorption at lower energy in the one sample 
that has been reported leaves its conduction state uncer- 
tain. The conducting state of CeN is uncertain. There 
are as yet no clear data for the rest of the series. 



3. Theoretical advances, band gaps 

We have already to some extent discussed the band 
gaps in the previous sections. Here we resummarize the 
results with this focus. 

The SIC calculations of Aerts et al.^ predicted both 
half-metallic (PrN-GdN), insulating (TbN-HoN) and 
semimctallic (ErN-YbN) behavior would occur in the 
REN series. As already mentioned, most of the LSDA+/7 
calculations find that when only Uf is included, the RENs 
have a small band overlap and are thus semimetals. This 
continues the trend of t he oth er RE-V, in particular the 
Gd pnictides and ErAsP^EH Base d on the early exper- 
imental workpID Larson et alP^ pointed out the need 
for a Ud parameter or shift of the d bands to open a 
gap. This seemed also justified based on Lambrecht's 
analysis^ of scaling the shifts inversely proportional to 
the dielectric constant from known Fermi surface prop- 
erties of the semimetallic pnictides. Using this approach, 
Larson et alp3 found a systematic upward trend of the 
direct band gap at X from LaN to LuN both in the ex- 
perimental data and in the calculations. This correlated 
with the decreasing lattice constants. However, the fluc- 
tuations of the individual values around this general be- 
havior are fairly large and it is not clear if these are due 
to experimental uncertainties or whether the LSDA+/7 
level of theory is able to capture these details. There is 
certainly no good agreement in the fluctuations of theory 
and experiment around the general increasing trend. One 
of the main successful predictions of the theories includ- 
ing / electrons explicitly and in contrast to the earlier 
models treating / electrons as core electrons, is the pre- 
diction of a gap red-shift in GdN, similar to that occur- 
ing in EuO. The clear experimental confirmation of this 
redshift came from Trodahl et aXP^ and led to a refine- 
ment of the Ud parameter for GdN. It is not clear however 
how good the approximation is to keep this Ud parameter 
constant across the REN series. The commonly made ap- 
proximation of averaging the majority and minority spin 
gap to obtain the effective gap above Tq was explicitly 
tested with non-collinear spin calculations with randomly 



pointing spins by Mitra. These results were in clud ed in 
Ref. M The HSE 120 and QSGW calculations^! both 
predict a very small gap in GdN when applied at the 
experimental lattice constant. 

The general band structure picture is thus well ac- 
cepted. An indirect gap very close to zero between 
majority-spin valence bands at T and majority-spin con- 
duction bands at X is followed by a larger minority-spin 
gap. In some calculations this gap is slightly negative 
and in others slightly positive depending on the details 
of the method. The optical absorption however is domi- 
nated by the smallest direct gap at X. Above the mag- 
netic ordering temperature, one presumably observes an 
average of the smaller majority-spin and higher minority- 
spin gaps. In earlier work that absorption threshold 
was placed at 0.98 eV in GdN, while in the more re- 
cent measurements it is found to be at about 1.3 eV. 
Below Tc, the spin-polarization results in a smaller gap 
between the majority-spin states. The theory and exper- 
iment agree pretty well on this redshift which is found 
to be 0.3-0.4 eV. These predictions are common to the 
LSDA+f// + U d theory, the HSE and the GW. For the 
other RENs, there are fewer detailed studies available to 
check experiment and theory. However, the anomalous 
band structure predicted by LSDA+J7 for EuN with a 
band crossing the Fermi level is not supported by the 
experiment or by the GW calculations. 

It would be quite interesting to find direct confirma- 
tion in the optical spectra of these different spin-up and 
spin-down bands but this has not yet been possible ex- 
perimentally. 



4- Theoretical advances, interband optical response 
functions 

Higher energy interband optical transitions have been 
measured for only YbNp31 but a theoretical treatment 
calculating the optical resp onse functions for the RENs 
has been reported by MitraP^l These could provide de- 
tailed information on the band structure when spectro- 
scopic ellipsometry data become available in the future. 
In particular, they analyse group theoretical selection 
rules and provide the only fully symmetry-labelled band 
structure for GdN. 

A study of optical response functions was also carried 
out by Ghosh et alp23 The two papers however differ 
rather strongly in details. Ghosh used two slightly dif- 
ferent LSDA+C7 implementations. One simply adds Uf. 
The other includes a smaller Uf but additional shifts of 
the d and / bands adjusted separately the position of the 
empty and filled / states with respect to the d-bands. 
This latter calculation gives a small gap in GdN but no 
attempt was made to adjust Ud to reproduce the exper- 
imental band gap. They report calculated optical con- 
ductivity, dielectric functions and reflectity for the entire 
series of Gd pnictides. For GdN, their optical conduc- 
tivity shows a broad spectrum peaking at about 8 eV, 
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similar to the results of Mitra and Lambrecht p^l which 
however, peak at 10 eV. They identify a doublet of peaks 
at 3-4 eV with N-p to Gd-d transitions and a peak at 
7 eV with occupied Gd-/ to empty Gd-d transitions. In 
contrast, Mitra and Lambrecht provide a much more de- 
tailed discussion of the fine structure of the dielectric 
function. They instead assign a peak at 7 eV in the di- 
electric function to a transition from the valence band to 
unoccupied / bands because they only see this peak in 
the minority-spin channel. On the other hand, they find 
no strong optical transitions from the occupied / bands 
to the conduction band. 

The refle ctivi ty spectra look very different, because in 
Ghosh et alP^ these are dominated by the strong peak 
at low energy due to the zero gap in their calculation. 
Although not very clearly stated, it appears that they 
include the intraband Drude contribution in these spec- 
tra which leads to a high reflectivity of nearly 100 %, 
dropping down to about 25 % within less than 1 eV. In 
contrast, Mitra's results, which do not include intraband 
transitions and has a finite gap, show a much lower re- 
flectivity at zero energy, also near 25%, but then shows 
the details of this peak at higher energies. In Ghosh's 
figures, these details are not visible in their results. The 
results of the reflectivity of GdAs and GdSb by Ghosh 
et alP^Sl agree well with experiment. Given that in real 
GdN there is always a large free carrier concentration, 
their calculated spectra may still be quite relevant for 
the overall expected behavior of the reflectivity. 



5. Magneto- optics 

Interestingly, Ghosh et alP^ also discuss the magneto- 
optic Kerr effect in the Gd compounds. Remark able 
magneto -opt ic results were earlier discovered in CeStP^ 
and Ce£p2! anc j successfully addr esse d with LSDA+C/ 
calculations by Liechtenstein et al.J^Obut remain to be 
studied experimentally in the nitrides. T he o ptical prop- 
erties of CeN were studied by Delin et alP^ and provide 
evidence for the 4/ band formation in CeN within the 
LSDA+f/ method. 



B. X-ray spectroscopies 

1. Experimental results 

In view of the considerable uncertainties of the theory, 
and the need to accurately locate the 4/ electron ener- 
gies, validation of the higher-energy electronic structure 
by experiment is important. Here again the majority of 
available X-ray photoelectron (XPS), absorption (XAS) 
and emission (XES) spectroscopic data are for GdN. XPS 
results by the Gottingen group provided the first infor- 
mation about the filled states, locating the Gd-4/ lev- 
els about 8 eV below the Fermi levelP^ The same study 
investigated the N-p projected density of empty states 



by XAS at the N K-edge, and the Gd-4/ empty states 
by Gd M-edge absorption. The XMCD signals at both 
edges were also investigated, as discussed above. The 
K-edge data were compared with computed results by 
Aerts et al.J^ providing the very first theory-experiment 
confrontation for the higher-energy structure in the con- 
duction band density of states, with the results serving 
to drive further experimental and theoretical work in the 
search for close agreement between the two. The M-edge 
data reflect primarily the atomic multiplet structure in 
the 4/ shell, showing no discernable solid-state effects. 
T he Gd £2,3 spectra were measured by Leuenberger et 
a ^| 59 | i44 | an( j com p arec i with t heor y by Abdelouahed and 
AlouanP^ and Antonov et alPSl 

Preston et al. repo rted XAS/XES spectra from epi- 
taxial GdN films The spectra were compared with 
N-p like partial densities of states from LSDA+J7 calcu- 
lations tuned to the measured optical gap by Larson et 
al.pSl obtaining reasonable agreement between the main 
features in both the filled and empty levels. The effects 
of core-holes were also investigated. The Gd M^-edges 
were also studied in both absorption and emission. This 
combination provided the first measurement of the 4/ 
filled-to-empty level splitting as about 12.5 eV, which, in 
conjunction with the XPS measurement of the occupied 
4/ levels' location, places the empty 4/ levels at about 
5 eV above the conduction band edge. 

More detailed XAS/XES spectra for SmN and DyN 
were reported by Preston et alP^ and compared to 
LSDA+f/ densities of states. The valence-to-conduction 
band gap was suggested by the results to lie close to 
1.5 and 1 eV for DyN and SmN, respectively, though 
the estimate is problematic due to the uncertain core- 
hole influence on the absolute energy scale of the XAS 
spectra. A full set of XPS, XAS, and XES spectra from 
HoN have b een recently reported to investigate p — / 
hybridi satio n^-^ Further XAS data have been reported 
for LuNliS and EuNp3 Figure 14 provides a compari- 
son between the N K-edge XAS d ata from the set of 
six RENs, along with metallic HfNP^ Band structure 
calculations and symmetry considerations show that the 
largest features, labelled A and B, correspond to crys- 
tal field split Gd-5d t^g and e g levels, respectively. This 
assignment is supported by the trend across the series, 
where the e g peak shifts upwards in energy as the lattice 
constant decreases reflecting the fact that these orbitals 
are orie nted directly towards the nearest neighbour N 
ions? 
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C. Phonon spectra 

The phonon properties of the RENs have also received 
some attention from theorists. In the rocksalt structure, 
there is no first order allowed Raman spectrum because 
of the presence of an inver sion cent er. Nonetheless, Ra- 
man spectra were 

observed™^ They were shown to 
be disorder induced Raman spectra and dominated by 
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FIG. 14. Nitrogen K-edge X-ray absorption spectra from a 
series of REN films (from Refs. \49\ ELM QSIl and [173} along 
with metallic HfN (Ref. ITgTj) . The peaks near 401 eV and 
405 eV are associated with crystal field split RE 5d t2 a and 
e 9 states, respectively. Spectra offset for clarity. 

the LO-phonons at the Brillouin zone L points. This 
vibrational mode corresponds to a breathing mode and 
is thus most strongly excited with above band gap reso- 
nant Raman excitation. The role of the various phonons 
was first elucidated in ScN, a transition metal cousin 
of the RENs by means of full phonon den sity of states 
and phonon band structure calculations!^ 1 In the RENs, 
pseudopotential calculations were not found to be suc- 
cessful and hence a frozen phonon method using the 
LSDA+f/ band structure approach^ for the total en- 
ergy was used. Phonons were calculated at T, X and 
L and the LO-phonons at L were found to account well 
for the observed Raman peaks. Their trends with lattice 
constant were explained and agree with the observed be- 
havior of lattice c onst ant as a function of the RE atomic 
number. Jha et alP^ studied the phonon dispersion and 
density of states in Yb pnictides, including YbN using a 
semi-empirical three-body potential. 

VII. SUMMARY 

The RENs were first investigated 50 years ago, but 
with results that were recognised to be subjected to 
both Vjv and reaction with oxygen. However, recent ad- 
vances in UHV-based thin-film growth and passivation 
have permitted substantially enhanced reliability of the 
data concerning their electronic and magnetic proper- 
ties, with results that suggest devices for both spintron- 
ics and as conventional narrow-gap semiconductors. The 
experimental work is still very incomplete; only a few of 
the series have even reasonably well-established conduct- 
ing states, semiconducting or metallic. However, even 
within those few there are some clear opportunities to 
form heterojunctions that present potentially exploitable 
properties; indeed there is already demonstrated a GdN- 



based spin filter.^ In concert with experimental inves- 
tigation these compounds have been studied even more 
thoroughly within various theoretical treatments. In that 
context they are attractive as a test for treatments of 
strong correlations in the presence of a relatively simple 
NaCl crystal structure. Thus much of this review sum- 
marises exactly the various treatments used in electronic 
and magnetic structure calculations, outlining the rela- 
tive strengths and weaknesses in both their fundamental 
underpinning and their agreement with experimental re- 
sults. The points of contact between theory and experi- 
ment lie primarily in the measured Curie temperatures, 
the magnitude of optical band gaps and the delineation 
of higher-energy excitation with X-ray absorption and 
emission spectroscopies. 

Recent thin film growth has been reported by physical 
vacuum deposition (CeN, SmN, EuN, GdN, DyN, HoN, 
ErN, LuN), sputtering (GdN) and CVD (GdN, DyN). 
In most cases the films have been subjected to mag- 
netic measurements, and for a smaller selection simple 
transport studies have indicated heavily doped ferromag- 
netic semiconductor behaviour. The most fundamental 
magnetic parameters (Curie temperature, saturation mo- 
ment) are in substantial agreement with the historical 
data, and are now supported by XMCD results for at 
least a small selection of the series. That agreement is 
in stark contrast to their electronic behaviour, with most 
so far investigated (SmN, GdN, DyN, ErN, and possi- 
bly EuN, HoN) showing clearly semiconductor-like tem- 
perature dependences of resistivity. Optical interband 
absorption and weak free-carrier absorption corroborates 
that conclusion in at least some of the REN series. 

Epitaxial film growth has so far been demonstrated 
for only three of the series; GdN, EuN and SmN. Even 
with this very limited set the recent data show a striking 
contrast; GdN has a huge magnetic moment and a very 
weak coercive field whereas SmN is a near-zero moment 
ferromagnet with a coercive field some 2-3 orders of mag- 
nitude larger. Coupled with their similar semiconductor 
electronic structure the combination of the two suggests 
potential as a magnetoresistance memory element. 

The RE elements, including its end members (La, Lu) , 
comprise 15 species. One, Pm, has no stable isotopes, 
making it inappropriate for any condensed matter inves- 
tigation. Nonetheless, there has been no more than a 
very dilute start at investigating the materials within 
the recent thin-film decade. Even for the most thor- 
oughly studied, GdN, there remains a vexing theoretical- 
experimental discrepancy concerning its Curie tempera- 
ture. There is a clear need for continued work to settle 
these and many other questions. 
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